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In Paper I of this series we described observations of 15 extragalactic targets taken 
with the Hubble Space Telescope+GHRS+G160M grating for studies of the low-z Lya 
forest. We reported the detection of 111 Lya absorbers at significance level (SL) > 3a, 
81 with SL > 4cr, in the redshift range 0.002 < z < 0.069, over a total pathlength cAz = 



(N 
> 

OO i 116, 000 km s~ x (Az = 0.387). In this second paper, we evaluate the physical properties 

of these Lya absorbers and compare them to their high-z counterparts. The SL > 4a 
distribution of Doppler parameters is similar to that at high redshift, with (b) = 38.0 ± 
15.7 km s _1 . The true Doppler parameter may be somewhat lower, owing to component 
blends and non-thermal velocities. The distribution of equivalent widths exhibits a 
significant break at W < 133 mA, with an increasing number of weak absorbers (10 mA 
< W < 100 mA). Adopting a curve of growth with b = 25 ± 5 km s 1 and applying 
a sensitivity correction as a function of equivalent width and wavelength, we derive 
the distribution in column density, jVHi~ L72±0 - 06 for 12.5 < log [N m ] < 14.0. We 
find no redshift evolution within the current sample at z < 0.07, but we do see a 
significant decline in d/V /dz compared to values at z > 1.6. Similiar to the high 
^ ■ equivalent width (W > 240 mA) absorbers, the number density of low-W absorbers at 

z=0 is well above the extrapolation of cLV /dz from z > 2, but we observe no difference 
in the mean evolution of d/V /dz between absorbers of high (W > 240 mA) and low 
(W < 100 mA) equivalent width absorbers. While previous work has suggested slower 
evolution in number density of lower- W absorbers, our new data do not support this 
conclusion. A consistent evolutionary pattern is that the slowing in the evolution of 
the low column density clouds occurs at lower redshift than for the higher column 
density clouds. A 4 — 5a signal in the two-point correlation function of Lya absorbers 
for velocity separations Av < 150 km s _1 is consistent with results at high-z, but 
with somewhat greater amplitude. Applying a photoionization correction, we find that 
the low-z Lya forest may contain ~ 20% of the total number of baryons, with closure 
parameter £liya = (O.OOSiO.OOl)/!^ 1 , for a standard absorber size and ionizing radiation 
field. Some of these clouds appear to be primordial matter, owing to the lack of detected 
metals (Si III) in a composite spectrum, although current limits on composite metallicity 
are not strong. 
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Subject headings: intergalactic medium — quasars: absorption lines — ultraviolet: 
galaxies 



1. Introduction 

Since the discovery of the high-redshift Lya forest over 25 years ago, these abundant absorption 
features in the spectra of QSOs have been used as evolutionary probes of the intergalactic medium 
(IGM), galactic halos, large-scale structure, and chemical evolution. Absorption in the Lya forest 
of H I (and He II) is also an important tool for studying the high-redshift universe (Miralda-Escude 
& Ostriker 1990; Shapiro, Giroux, & Babul 1994; Fardal, Giroux, & Shull 1998). A comparison of 
the H I and He II absorption lines provides constraints on the photoionizing background radiation, 
on the history of structure formation, and on internal conditions in the Lya clouds. In the past 
few years, these discrete Lya lines have been interpreted in the context of N-body hydro dynamical 
models (Cen et al. 1994; Hernquist et al. 1996; Zhang et al. 1997; Dave et al. 1999) as arising from 
baryon density fluctuations associated with gravitational instability during structure formation. 
The effects of hydrodynamic shocks, Hubble expansion, photoelectric heating by AGN, and galactic 
outflows and metal enrichment from early star formation must all be considered in understanding 
the IGM. 

At high redshift, the rapid evolution in the distribution of Lya absorption lines per unit red- 
shift, dj\f jdz oc (1 + z) 7 (7 f» 2.5 for z > 1.6), was consistent with a picture of these features as 
highly ionized "clouds" whose numbers and sizes were controlled by the evolution of the IGM pres- 
sure, the metagalactic ionizing radiation field, and galaxy formation. However, the rapid evolution 
of the Lya forest stopped at z < 1.6. One of the delightful spectroscopic surprises from the Hubble 
Space Telescope (HST) was the discovery of Lya absorption lines toward the quasar 3C 273 at 
z cm = 0.158 by both the Faint Object Spectrograph (FOS, Bahcall et al. 1991) and the Goddard 
High Resolution Spectrograph (GHRS, Morris et al. 1991, 1993). The number of these absorbers 
was found to be far in excess of their expected number based upon an extrapolation from high-z 
(see e.g. Weymann et al. (1998) and section 5 therein). This evolutionary shift is probably a result 
of the collapse and assembly of baryonic structures in the IGM (Dave et al. 1999) together with the 
decline in the intensity of the ionizing radiation field (Haardt & Madau 1996; Shull et al. 1999b). 
Detailed results of the Lya forest evolution in the redshift interval < z < 1.5 are described in the 
FOS Key Project papers: the three catalog papers (Bahcall et al. 1993, 1996; Jannuzi et al. 1998) 
and the evolutionary analysis study (Weymann et al. 1998). 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract 
No. NAS5-26555. 

2 Also at JILA, University of Colorado and National Institute of Standards and Technology. 
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However, the HST/FOS studies were primarily probes of strong Lya lines with equivalent 
widths greater than 0.24 A. A great deal more information about the low-z IGM can be gained 
from studies of the more plentiful weak absorbers. Realizing the importance of spectral resolution 
in detecting weak Lya absorbers, the Colorado group has engaged in a long-term program with 
the HST/GHRS, using the G160M grating at 19 km s -1 resolution, to study the very low-redshift 
(z < 0.07) Lya forest. Earlier results from our study have appeared in various research papers 
(Stocke et al. 1995; Shull, Stocke, & Penton 1996; Shull et al. 1998) and reviews (Shull 1997; 
Shull, Penton, & Stocke 1999a). 

The current series of papers discusses our full GHRS dataset. In Paper I (Penton, Shull, & 
Stocke 2000a) we described our HST/GHRS observations and catalog of Lya absorbers. In Paper 
II (this article) we describe the physical results from our program, including information on the 
physical parameters and nature of the low-redshift Lya forest. A discussion of the relationship 
between the Lya absorbers discovered in our GHRS program and the nearby galaxy distribution as 
mapped using available galaxy redshift survey data will be presented in Paper III (Penton, Stocke, 
& Shull 2000b). We believe that the low-redshift Lya forest of absorption lines, combined with 
information about the distribution of nearest galaxies, can provide a probe of large-scale baryonic 
structures in the IGM, some of which may be remnants of physical conditions set up during galaxy 
formation. In Tables 1 and 2 we present the basic data for the definite (> 4cr) and possible 
(3 — 4cr ) Lya absorbers from Paper I. In general, scientific results will be determined for only the 
> 4a absorbers in Table 1, with data using the expanded (> 3a) sample (Tables 1 and 2) shown 
only as corroborating. Only the absorbers judged to be intervening ( "intergalactic" ) are included 
; see Paper I for detailed criteria and analysis. We exclude all Galactic metal-line absorbers and 
absorbers "intrinsic" to the AGN, with (czagn — c z abs) < 1200 km s _1 ; see Paper I. The information 
in Tables 1 and 2 by column is: (1) name of target AGN; (2) absorber wavelength and error in A; 
(3) absorber recession velocity and error quoted non-relativistically as cz in km s _1 ; (4) observed b- 
value (fe bs) and error in km s _1 ; (5) resolution-corrected b- value and error in km s" 1 ; (6) rest-frame 
equivalent width and error in mA; (7-10) estimated column densities in cm -2 assuming 6-values of 
20, 25, 30 km s _1 and the value from column (5). Detailed descriptions of the determination of 
values in columns (1-6) can be found in Paper I. As described in detail in Paper I, the uncertainties 
for the W values in column 6 are the uncertainties in the Gaussian fit to each feature and not 
the significance level (i.e., typically W/AW > significance level (SL)). Further discussion of the 
correction of the 6 b s -values (columns 4 and 5) can be found in § 2. Descriptions and justification 
for values in columns (7-10) can be found in § 3. 

In this paper, we analyze physical quantities derivable from the measured properties of the 
intergalactic low-z Lya lines of Paper I. In § 2, we discuss the results and limitations of the b- 
value determinations for our low-z Lya detections. In § 3, we discuss the basic properties of our 
measured rest-frame equivalent width (W) distributions and compare them to higher- z distribu- 
tions. In § 4, we estimate H I column densities (Ahi) for our Lya absorbers and discuss their 
distribution, (&Sf /dz) z= o, relative to similarly derived values at higher redshift. In § 5, we discuss 
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the z distribution of the low- z Lya forest, as well as the cumulative Lyman continuum opacity of 
these absorbers and the z evolution of the number density of lines, <\M /dz. In § 6, we analyze the 
cloud-cloud two-point correlation function (TPCF) for low-z Lya clouds, and in § 7, we explore 
the limits on metallicity of the low-z Lya forest. In § 8, we estimate its contribution, Vtiya-, to the 
closure parameter of the Universe in baryons, Ob, inferred from D/H. Section 9 summarizes the 
important conclusions of this investigation. 



2. Observed b- value Distribution 

Doppler widths (6-values) are estimated from the velocity widths (Wg= b h s /^/2) of our Gaus- 
sian component fits. As such, they are not true measurements of the actual 6-values, as when 
fitting Voigt profiles, but rather velocity dispersions assuming that the absorption lines are not 
heavily saturated. This is a particularly good assumption for the large number of low-W lines (i.e., 
W < 75 mA), but it becomes increasingly suspect for the higher W lines. 

The GHRS/G160M produced spectral resolution elements (REs) with Full Widths at Half 
Maximum (FWHM) of ~19 km s _1 . The line spread function (LSF) of the HST+GHRS/G160M 
is approximately Gaussian with oq ~ 8.06 km s^ 1 for both pre- and post-COSTAR data (Gilliland 
et al. 1992; Gilliland & Hulbert 1993). As discussed in Paper I, to improve the robustness of 
our Gaussian component fitting we smooth our data with the LSF. The measured b- value (fe bs) 
is actually the convolution of the instrumental profile, our pre-fit smoothing, and the observed 
b- value (6) of the absorber. The 6-values add in quadrature, b 2 = bl bs — 2 * (Tq = 2(Wg — <Tq), 
where Wg is the Gaussian width of the fitted absorption component. Therefore, we are hampered 
in detecting absorptions with 6-values near or below b = \[2 oq = 11.4 km s _1 . Tables 1 and 2 
present 6- values, rest frame equivalent widths (W), and estimated H I column densities for our 
SL > 4a and 3a < SL < 4a Lya samples, respectively. The H I column densities are estimated 
assuming b = 20, 25, 30 km s -1 , and the individual corrected b- value for each absorber. 

Motion of the target in the GHRS large science aperture (LSA) during our lengthy exposures 
can broaden the line spread function (LSF) or modify the wavelength scale of our spectra, causing 
us to overestimate the 6-values. Our subexposures are generally of insufficient signal-to-noise to 
perform accurate cross-correlations to minimize the wavelength scale degradation. Weymann et al. 
(1995) performed such cross-correlations on one GHRS G160M spectrum of 3C 273 in our sample. 
For the 1220 A and 1222 A features in 3C 273, they measured 5-values of 40.7 ± 3.0 km s" 1 
and 34.3 ± 3.3 km s _1 , respectively. As indicated in the first two entries of Table 1, we measure 
much larger b- values of 69±5 and 72±4 km s -1 . We believe these differences arise from target 
motion within the LSA, causing larger 6-values, although the line center and W measurements are 
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Table 1: H I Column Densities (Nm) of Definite Lya features 
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13.52 
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Note. — * Indicates that this absorber was observed with HST/GHRS pre-COSTAR. 
a b- value after correcting for the GHRS intrumental profile and our pre-fit smoothing (see § 2). 
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Table 2: H I Column Densities (Nui) of Possible Lya features 



Target 


Wavelength 


Velocity 


bobs 




VV 




log[JV H i 


(in cm 


)] 




A 




km s~ 1 


km s 1 


km s 


-l 


mA 


6=20 


b=25 


b=30 


fa— col 5 a 


.30273 


1224.587 


± 


0.150 


2199 ± 37 


59 ± 53 


57 ± 


55 


29 ± 35 


12.77 


12.76 


12.75 


12.74 


*3C273 


1234.704 


± 


0.029 


4694 ± 7 


71 ± 60 


69 ± 


61 


25 ± 29 


12.69 


12.68 


12.68 


12.67 


*3C273 


1265.701 


± 


0.064 


12338 ± 16 


37 ± 24 


33 ± 


26 


21 ± 18 


12.61 


12.61 


12.60 


12.60 


*3C273 


1266.724 


± 


0.084 


12590 ± 21 


54 ± 29 


52 ± 


30 


24 ± 18 


12.68 


12.67 


12.67 


12.66 


*3C273 


1268.969 


± 


0.076 


13144 ± 19 


45 ± 26 


43 ± 


28 


18 ± 15 


12.55 


12.55 


12.55 


12.54 


AKN120 


1223.088 


± 


0.039 


1829 ± 10 


30 ± 14 


25 ± 


17 


64 ± 48 


13.15 


13.13 


13.12 


13.13 


AKN120 


1247.267 


± 




7792 ± 26 


38 ± 29 


34 ± 


32 


19 ± 22 


12.57 


12.57 


12.56 


12.56 


ES0141-G55 


1249.932 


± 


0.036 


8449 ± 9 


23 ± 11 


17 ± 


15 


12 ± 8 


12.36 


12.36 


12.36 


12.37 


ES0141-G55 


1252.483 


± 


0.041 


9078 ± 10 


28 ± 13 


23 ± 


16 


12 ± 7 


12.35 


12.34 


12.34 


12.34 


FAIRALL9 


1265.407 


± 


0.029 


12265 ± 7 


25 ± 13 


19 ± 


18 


11 ± 8 


12.31 


12.31 


12.31 


12.31 


H1S21+643 


1238.014 


± 


0.036 


5510 ± 9 


26 ± 11 


21 ± 


14 


23 ± 13 


12.66 


12.65 


12.65 


12.66 


H1821+643 


1240.569 


± 


0.036 


6140 ± 9 


21 ± 10 


14 ± 


16 


24 ± 16 


12.66 


12.66 


12.65 


12.67 


H1S21+643 


1244.966 


± 


0.031 


7225 ± 8 


22 ± 8 


15 ± 


12 


25 ± 13 


12.69 


12.68 


12.68 


12.70 


H1S21+643 


1247.362 


± 


0.029 


7815 ± 7 


39 ± 31 


36 ± 


34 


27 ± 29 


12.73 


12.73 


12.72 


12.72 


H1S21+643 


1252.477 


± 


0.042 


9077 ± 10 


26 ± 12 


21 ± 


15 


23 ± 15 


12.65 


12.65 


12.64 


12.65 


H1S21+643 


1254.874 


± 


0.099 


9668 ± 24 


29 ± 26 


24 ± 


31 


21 ± 25 


12.61 


12.61 


12.61 


12.61 


MARK279 


1237.915 


± 


0.029 


5486 ± 7 


34 ± 26 


30 ± 


29 


17 ± 18 


12.53 


12.52 


12.52 


12.52 


MARK279 


1238.502 


± 


0.047 


5631 ± 12 


26 ± 19 


21 ± 


23 


18 ± 18 


12.54 


12.53 


12.53 


12.54 


MARK290 


1232.797 


± 


0.064 


4224 ± 16 


27 ± 23 


21 ± 


29 


41 ± 49 


12.93 


12.92 


12.91 


12.93 


MARK290 


1235.764 


± 


0.044 


4955 ± 11 


30 ± 15 


26 ± 


18 


28 ± 19 


12.74 


12.73 


12.73 


12.73 


MARK290 


1245.869 


± 


0.026 


7447 ± 6 


18 ± 5 


8 ± 


12 


18 ± 7 


12.54 


12.54 


12.53 


12.55 


*PKS2155-304 


1234.767 


± 


0.051 


4709 ± 12 


36 ± 21 


32 ± 24 


15 ± 14 


12.45 


12.45 


12.44 


12.44 


*PKS2155-304 


1246.990 


± 


0.029 


7724 ± 7 


35 ± 31 


31 ± 35 


13 ± 16 


12.40 


12.39 


12.39 


12.39 


*PKS2155-304 


1247.510 


± 


0.029 


7852 ± 7 


34 ± 27 


30 ± 31 


13 ± 16 


12.41 


12.40 


12.40 


12.40 


*PKS2155-304 


1255.084 


± 


0.041 


9720 ± 10 


29 ± 13 


24 ± 16 


13 ± 9 


12.41 


12.40 


12.40 


12.40 


*PKS2155-304 


1256.636 


± 


0.042 


10102 ± 10 


30 ± 13 


25 ± 15 


14 ± 8 


12.42 


12.41 


12.41 


12.41 


PKS2155-304 


1264.806 


± 


0.058 


12117 ± 14 


42 ± 19 


39 ± 20 


31 ± 19 


12.78 


12.78 


12.77 


12.77 


Q1230+0115 


1236.045 


± 


0.041 


5025 ± 10 


29 ± 12 


24 ± 15 


53 ± 31 


13.05 


13.04 


13.03 


13.04 


Q1230+0115 


1242.897 


± 


0.044 


6714 ± 11 


25 ± 14 


19 ± 19 


45 ± 35 


12.97 


12.96 


12.95 


12.98 


Q1230+0115 


1246.254 


± 


0.049 


7542 ±12 


35 ± 16 


31 ± 18 


44 ± 27 


12.96 


12.95 


12.95 


12.94 



Note. — * Indicates that this absorber was observed with HST/GHRS pre-COSTAR. 
a b- value after correcting for the GHRS intrumental profile and our pre- fit smoothing (see § 2). 

unaffected. We use the Weymann et al. (1995) 6-values for these two features to compute A^hi in 
column 10 of Tables 1 and 2, instead of our values, although these are pre-COSTAR values which 
will not be used statistically herein. Other spectra, such as those of Markarian 335, Markarian 501, 
and PKS 2155-304 (pre-COSTAR only), also seem to be affected by this degradation ("jitter"). 
The fact that all of the suspected exposures were taken before the 1993 HST servicing mission 
leads us to speculate that spacecraft jitter, due to wobble introduced by thermal gradients across 
the solar panels during passage across the terminator, may be responsible for the spectral smearing 
we observe. In addition to installing COSTAR, this servicing mission corrected the jitter problem 
caused by spacecraft wobble (Bely, Lupie, & Hershey 1993; Brown 1993). Motions in the target 
aperture can cause spectral motion on the detector that can broaden spectral features and increase 
measured b- values (an offset of 1 arcsec in the LSA corresponds to ~70 km s" 1 on the spectrum). 
All wavelengths in Tables 1 and 2 are LSR values, and all velocities are non-relativistic values 
relative to the Galactic LSR, as explained in Paper I. 

Figure 1 shows our 6-value distribution. Grey boxes indicate definite (SL > 4a) Lya 6-values, 
while the black boxes indicate possible (3a < SL < 4a) detections. Together, these two samples 
form our "expanded" sample (SL > 3a). The b- value distribution in the top panel of Figure 1 
shows a secondary peak at 47.5 < b < 72.5 km s _1 not present in the other observed 6-value 
distributions. This peak is produced primarily (16 of 22 absorbers, or 73%) by absorption features 
in pre-COSTAR data that appear to suffer from systematic broadening. In addition, we interpret a 
second population (5 of 22 absorbers, or 23%) as possible blended absorption features that we were 
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not able to resolve fully. One post-COSTAR absorption feature appears unblended and thus truly- 
broad. These blended features will be discussed further in § 6 when we analyze the linear two-point 
correlation of our absorbers. The second panel of Figure 1 shows only those 6-values obtained after 
the 1993 HST servicing mission. 

In Table 3, we report statistics for 6, with and without the pre-COSTAR absorbers included in 
the sample. Table 3 contains the following information by column: (1) the reference papers for the 
6-value sample; (2) the redshift range over which 6-values were determined; (3) the mean redshift 
of the absorber sample; (4) the observed wavelength range in A from which the absorber data were 
extracted; (5) the spectral resolution of the observations in km s _1 ; (6) the median 6-value of the 
sample in km s _1 ; (7) the mean 6-value of the sample in km s _1 ; and (8) the standard deviation 
of the mean in km s _1 . Where the referenced work reports two samples, the 6-value results in 
columns (6-8) are reported for both samples (one in parentheses), with the two samples described 
in footnotes to Table 3. Our post-COSTAR definite (SL > 4a) distribution of 6 has a mean of 
38.0 ± 15.7 km s _1 and a median of 34.8 km s -1 , while our pre-COSTAR definite distribution has a 
mean of 58. 6± 15.9 km s _1 and a median of 60.6 km s _1 . Table 3 also compares our results with other 
studies at higher redshift and is ordered by decreasing mean redshift, (z). For studies up to (z)= 
3.7, Table 3 reports the redshift range, (z), wavelength (A) range, resolution element (RE), median 
6, mean 6, and standard deviation (<r&) of the mean 6 for the various Lya absorption samples. The 
third through seventh panels of Figure 1 show the higher-z (1.6 <z<4.0) 6 distributions from Lu 
et al. (1996), Carswell et al. (1984), Carswell et al. (1991), Kulkarni et al. (1996), and Khare et al. 
(1997). Statistics for these distributions are tabulated in Table 3. Lu et al. (1996) 's observations 
of the Lya-forest at 3.4 < z < 4.0, shown in the third panel of Figure 3, were taken with the Keck 
HIRES spectrograph, with a resolution element of 6.6 km s _1 . As expected, the Lu et al. (1996) 
distribution extends to lower 6-values than ours. This high-z distribution contains 412 absorbers 
with a peak at 6=23 km s _1 and a median value of 6=27.5 km s -1 . Like the high-z distribution 
of Lu et al. (1996), the low-z distributions show an increasing number of absorbers at decreasing 
6-values, until one approaches the resolution limit. Lu et al. (1996) report a detectable turnover 
in the 6-value distribution below 20 km s _1 with almost no absorbers at 6 < 15 km s -1 . Because 
these Keck data have three times higher resolution than HST/GHRS, features that are blended in 
our data would be resolved by Keck. 

A more accurate comparison to our data can be made with 6-value distributions taken at 
comparable resolution (e.g., Carswell et al. 1984). As indicated in Table 3, the 6-values of Carswell 
et al. (1984); Kulkarni et al. (1996); Khare et al. (1997) were obtained with spectral resolutions 
similar to the GHRS/G160M (~19 km s _1 ). For direct comparison to our 6-values, the second- 
to-last panel in Figure 1 displays the combined 6-value distributions of these three studies and is 
labeled "C84KK" . The final panel of Figure 1 displays the combined 6-value distributions of other 
studies (Lu et al. 1996; Carswell et al. 1984, 1991; Kulkarni et al. 1996; Khare et al. 1997) and is 
labeled "All" . Comparing our 6-value distributions to the samples listed here, we conclude that at 
the 2a level, the high and low-z 6-value distributions appear to differ, with the low-z distribution 
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Table 3. Comparison of our b- values to other Lya studies 



Reference 


z range 




A range 


Resolution 


Median b 


Mean b 


_ a 
Ob 








(A) 


( km s" 1 ) 


( km s" 1 ) 


( km s" 1 ) 


( km s" 1 ) 


Lu et al. (1996) b 


3.43-3.98 


3.70 


5380-6050 


6.6 


27.5(25.9) 


34.4(30.6) 


29.0(15.3) 


Kim et al. (1997) c 


3.20-3.51 


3.36 


5105-5484 


8 


30(27) 


NR d 


NR 


Hu et al. (1995) 


2.54-3.20 


2.87 


4300-5100 


8 


-35 


28 


10 


Rauch et al. (1992) 


2.32-3.40 


2.86 


4040-5350 


23 


33 


36 


16 


Kim et al. (1997) d 


2.71-3.00 


2.86 


4510-4863 


8 


35.5(30.0) 


NR 


NR 


Kirkman & Tytler (1997) 


2.45-3.05 


2.75 


4190-4925 


7.9 


~28 


23 


14 


Rauch et al. (1993) 


2.10-2.59 


2.34 


3760-4360 


~8 


26.4 


27.5 


NR 


Kim et al. (1997) d 


2.17-2.45 


2.31 


3850-4195 


8 


37.7(31.6) 


NR 


NR 


Carswell et al. (1984) 


1.87-2.65 


2.26 


3490-4440 


19 


25.0 


27.9 


10.9 


Khare et al. (1997)° 


1.57-2.70 


2.13 


3130-4500 


18 


27.7(30.6) 


29.4(32.1) 


7.9 


Carswell et al. (1991) 


1.84-2.15 


1.99 


3434-3906 


9 


33.0 


34.3 


14.1 


Kulkarni et al. (1996) f 


1.67-2.10 


1.88 


3246-3769 


18 


31.6(29.5) 


35.6(31.4) 


~15 


Savaglio et al. (1999) 


1.20-2.21 


1.70 


2670-3900 


8.5-50 g 


28.2 


31.1 


16.6) 


This Paper (SL > 4cr) h 


0.002-0.069 


0.035 


1218-1300 


19 


34.8(40.7) 


38.0(45.4) 


15.7(18.6) 


This Paper (SL > 3cr) h 


0.002-0.069 


0.035 


1218-1300 


19 


31.7(35.4) 


33.9(40.8) 


15.4(18.8) 



a The standard deviation of the mean, if available. 

b Lu et al. (1996) report values for two samples, A and B. The A sample contains the full set of detected Lya lines, 
while the B sample contains only those Lya lines whose values of N K j and 6 are well determined. Their results are 
reported as A(B). 

c Kim et al. (1997) report values for two jV H i ranges 13.8 <log[iV H i ] <16.0 and 13.1 <log[iV H i ] <14.0. Their results 
are reported with the second, lower TVhi range in parentheses. The Kim et al. (1997) sample was divided into three 
2-ranges, with the (z) — 2.31 and (z)— 2.86 ranges reporting the combined results of two sightlines each. The (z) — 
3.36 entry reports the observations of a single sightline. 

d NR = Not Reported 

e Khare et al. (1997) report values for two samples, A and B. The A sample contains the set of detected Lya lines 
that do not show unusual profiles. Their B sample includes the full set of detected Lya lines. Their results are 
reported as A(B). 

'Kulkarni et al. (1996) report values for two samples, A and B. The A sample contains a conservative number of 
sub-components in blended lines, while the B sample is more liberal in allowing sub-components. Their results are 
reported as A(B). 

s Savaglio et al. (1999) use a combination of HST/STIS high and medium resolution data, combined with 
UCLES/AAT ground-based observations. 

h We divide both our SL > 4cr and SL > 3a samples into two components. Our A sample does not include pre- 
COSTAR 6-values since we have found them to be problematic. Our B sample contains all b-values. Our results are 
reported as A(B). 
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containing somewhat broader lines. Comparisons via KS tests to determine if our SL > 4cr post- 
COSTAR distribution is drawn from the same parent sample as the other and combined samples 
range from highly unlikely (C84KK, 1% and All, 0.3%) to probable (Carswell 1991, 30%). However, 
due to the variety of methods used in the observations, data reduction, and profile fitting, we cannot 
be certain that the observed modest 6-value evolution is real at this time. A consistent dataset 
with higher spectral resolution and high S/N is required; HST/STIS echelle data would be ideal. 

It has been proposed (Kim et al. 1997) that there is evolution of 6-values with redshift. Al- 
though the three Kim et al. (1997) entries in Table 3 do show a tendency to higher 6- values (broader 
absorption features) at lower z, this is not confirmed by the other studies listed in the table. Our 
results, at significantly lower z than any other study in Table 3, do not support significant evolu- 
tion. Neither does the 6-value distribution found at 1.2 < z < 2.2 by Savaglio et al. (1999) in a 
single sightline show any signs of the evolution suggested by Kim et al. (1997). Instead, a small 
decrease in 6-values from (z)= 3 to (z}= is found by "effective equation of state" models of the 
IGM fitted to Lya line- width data (Ricotti, Gnedin, & Shull 2000; Bryan & Machacek 2000). 

The measurement of the 6-value distribution is different from the measurement of other distri- 
butions such as W and AT HI , because it is difficult to correct for incompleteness (what we cannot 
detect). For example, consider a portion of a spectrum with a 4cr detection limit of 100 mA. 
When calculating the W distribution, we cannot detect features in this region with W < 100 mA. 
Therefore, we eliminate this pathlength for W < 100 mA when determining the true rest-frame 
equivalent width (number-density) distribution, ftf(W)/Az(W). Here, Az(W) is the available red- 
shift pathlength for detection of absorption features as a function of W. But, in our hypothetical 
spectral region, we cannot distinguish between a non-detection of a narrow absorption feature of 
6=20 km s _1 and a broad one of 6=80 km s -1 , if both produce features of W < 100 mA. In 
addition, it is possible to misidentify very broad, low-W features as continuum undulations. 

The accurate measurement of 6-values is important in determining the actual H I column 
densities (Ahi) of the saturated Lya absorbers, since each 6-value produces a different curve of 
growth for the upper range of column densities in Tables 1 and 2. The only reliable method of 
deriving 6-values for such weak lines is from higher Lyman lines such as Ly/3 (Hurwitz et al. 1998; 
Shull et al. 2000). For example, Hurwitz et al. (1998) found that the Ly/3 absorption strength 
observed by ORFEUS II was in strong disagreement with the predicted values based on two Lya 
Voigt profiles in the GHRS spectrum of 3C 273 sightline associated with the intracluster gas of 
the Virgo supercluster (Morris et al. 1991; Weymann et al. 1995). Hurwitz et al. (1998) observed 
Ly/3 at 1029.11 and 1031.14 A, corresponding to the first two Lya features of Table 1, with Ly/3 
equivalent widths of 145 ± 37 mA and 241 ± 32 mA, respectively. Weymann et al. (1995) measured 
the Lya 6-values of these features as 40.7 ±3.0 km s _1 and 34.3 ±3.3 km s _1 , respectively. The Ly(3 
observations of Hurwitz et al. (1998) imply 6-values as low as 12 km s -1 . Similar trends towards 
lower 6-values are found in initial Ly/3 studies with the Far Ultraviolet Spectroscopic Explorer 
(FUSE) satellite (Shull et al. 2000). This disagreement in estimating 6-values can be understood 
if the Lya absorption profiles include non-thermal broadening from cosmological expansion and 
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infall, and multiple unresolved components. Hu et al. (1995) reach the same conclusion, based 
upon z~3 Keck HIRES QSO spectra. The analysis of Hu et al. (1995) implies that the average 
z~3 Lya absorber could be well represented by ~3 components, each having 6~15 km s _1 . 

Therefore, we are faced with a dilemma: should we use our approximate 6-values inferred from 
the line widths, a constant value based upon the better known high-z distribution, or an adjusted 
constant value considering the Ly/3 results ? In § 4 we compare the Nui results from these three 
alternatives and discuss the merits of each. 
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Fig. 1. — Top panel (PI) shows the low-z distribution of line widths, b. Two features in the 
3C 273 spectrum have been adjusted in accordance with Weymann et al. (1995). P2 shows only 
6-values obtained with HST+COSTAR. In PI and P2, grey and black boxes indicate our SL > 4cr 
and 3a < SL < 4cr samples, respectively. P3-P7 show higher-z distributions from various sources 
(Table 3). Panel P8,"C84KK", combines all the higher-z distributions with resolutions comparable 
to ours (~19 km s^ 1 ). P9, "All", combines panels P3-P7. 



-12- 



3. Observed Rest-Frame Equivalent Width Distribution 

In Figure 2 we display the rest-frame equivalent width (W) distribution for all of our detected 
Lya features. The solid grey boxes in Figure 2 represent our definite (SL > 4a) sample, while the 
black boxes represent our possible (3a < SL < 4a) sample. As expected, we detect an increasing 
number of absorbers at decreasing W, down to our detection limit. As discussed in the previous 
section, our spectra are of varying sensitivity and wavelength coverage. This observed W distri- 
bution is not the true W distribution, because we have not yet accounted for incompleteness. To 
determine the true W distribution, we must normalize the W density by the available pathlength 
Az(W). 

The pathlength, Az, is actually a function of W and z, since our spectra have varying 4a 
detection limits across the waveband and each spectrum covers a different waveband (z range). 
Without this important sensitivity correction, Az(W,z), any interpretation of the W distribution 
is premature. Previous studies, such as the HST/FOS Key project (Bahcall et al. 1993), avoided 
this problem by considering the W distribution only above a universal minimum W detectable in 
all portions of all spectra in the sample. This forces Az(W,z) to be a constant pathlength, so 
that the W distribution is a true representation of the detected absorbers. However, this procedure 
eliminates information about the lowest, most numerous, W absorbers. Because our GHRS program 
probes the lowest equivalent widths of any low-z Lya program to date, it is important to use all 
available data for the weakest absorbers. In this section, we will explore the W distribution without 
regard to z and neglect any evolution of W with z (we find no evolution in numbers of Lya absorbers 
with redshift over our small Az range). We will examine the z distribution (0.002 < z < 0.069) 
and the W evolution with z in later sections. 

In Paper I, we discussed the available pathlength for each sightline in detail. In Figure 3, we 
display the available pathlength (Az) in terms of W for each spectrum. The solid line in each panel 
indicates the full observational pathlength, uncorrected for the regions of spectra not available for 
Lya detection due to Galactic, HVC, intrinsic, and non-Lya intervening absorption lines as well as 
our "proximity effect" limit. This "proximity limit" (czagn — cz a h s > 1, 200 km s" 1 ) eliminates 
potential absorption systems "intrinsic" to the AGN (see Paper I for a detailed description and 
justification). The dot-dashed line indicates the "effective" or available pathlength after we remove 
the spectral regions unavailable for Lya detection. In Figure 3, cAz(W) is given in units of 
Megameters per second (Mm s" 1 or 10 3 km s _1 ). Note that three targets (I ZW 1, Q1230+0115, 
and Markarian 501) have significantly poorer sensitivity than the rest of our sample and should 
be re-observed with HST, using the Space Telescope Imaging Spectrograph (STIS) or the Cosmic 
Origins Spectrograph (COS). In Figure 4, we display the cumulative distribution of Az(W) for all 
spectra, for SL > 4a detections. In this cumulative plot, the solid line indicates the uncorrected 
pathlength. Also indicated is the pathlength available for Lya detection after correcting for the 
spectral obscuration due to Galactic and HVC features, intrinsic and non-Lya intervening absorbers 
(Intrinsic+Non-Lya), and our "proximity limit". Indicated by the lowest (dotted) line is the 
cumulative available pathlength, for SL > 4a detections, after the indicated corrections have been 
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Fig. 2. — Observed distribution in rest-frame equivalent width (W) of our detected Lya absorbers. 
The solid grey boxes represent definite (SL > 4a) detections, while the black boxes indicate possible 
(3a < SL < 4a) absorbers. The W bins are 12 mA in width. This distribution has not been 
corrected for the non-uniform wavelength and sensitivity coverage of our observations. No intrinsic 
Lya absorptions are included in this distribution. 

made. The left axis of Figure 4 indicates the available pathlength (cAz) in units of Mm s -1 , 
while the right axis indicates Az. Our maximum available pathlength (cAz) is 116 Mm s _1 (or 
Az=0.387) for all features with W > 150 mA. Also indicated in Figure 4 is our TVhi pathlength 
availability for a constant 6=25 km s _1 . This pathlength compares to Az ~ 8 at z < 0.3 for the 
Key Project for strong absorbers (W > 240 mA). 



3.1. The Low-z Equivalent Width Distribution 

In Figure 5, we apply the effective pathlength correction of Figure 4 to the jV(VV) distribution 
of Figure 2. This procedure gives the true detected W number density, 

n(Wi) = M{Wi)/Az{Wi)AW « d 2 M/dz dW \ m , (1) 

corrected for the pathlength, Az(Wj), available to detect features at each Wj, without regard to z. 
The approximation that n(W) is equal to dAf /dW is limited by our sample size (jV) and by our 
bin size (AW= 12 mA). However, this approximation is more than sufficient for our purposes. It 
has been shown (Sargent et al. 1980; Young, Sargent, & Boksenberg 1982; Murdoch et al. 1986; 
Weymann et al. 1998) that d 2 M/dz dW, the pathlength-corrected number density of lines with 
respect to W, is well modeled by n(W) = (Cw/W* ) exp (— IV/IV* ), where W* is the characteristic 
equivalent width and Cw is the characteristic line density per unit redshift. In Figure 6, n(W) is 
plotted versus W in natural logarithmic form, for both our definite (SL > 4a) and our expanded 
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Fig. 3. — Each panel shows the available pathlength, cAz(W) in Mm s _1 , for a sightline as a 
function of 4cr rest-frame equivalent width (W). The solid line is the full cAz, uncorrected for 
known Galactic, HVCs, intrinsic, or intervening features. Regions of FWHMx2 of each >3a non- 
Lya absorber are removed from the available cAz for detecting Lya absorbers. The dashed line 
is the corrected cAz, including the removal of 10 pixels (~0.2 A) at each spectral edge and the 
application of our cz cm - 1,200 km s _1 "proximity limit". 
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Fig. 4. — Cumulative available pathlength (Az) as a function of W, for SX > 4a detections. The 
left axis gives non-relativistic cAz(W) in units of Mm s _1 , while the right axis gives Az. The 
solid (upper) line indicates the observed, uncorrected, Az(W). Dashed, dot-dashed, and dot-dot- 
dashed lines indicate cAz(W) after correcting for spectral obscuration due to non-Lya and intrinsic 
absorption systems (Intrinsic+Non-Lya) , Galactic and HVC absorption systems (Galactic+HVC), 
and our cz cm - 1,200 km s _1 "proximity" limit, respectively and individually. The dotted line 
indicates the fully corrected "effective" Az(W), after the indicated spectral regions unavailable for 
intervening Lya detection have been removed. Top axis indicates the log [Nja] corresponding to 
W for b = 25 km s _1 ; for example, log [Nm] = 12.26 for W = 10 mA (left-hand border of bottom 
axis) . 
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Fig. 5. — Rest-frame equivalent width (W) number distribution, n(W), corrected for available 
pathlength (Az). The corrected number density, n(W)=J\f(W)/Az(W), is the observed W distri- 
bution corrected for the non-uniform sensitivity function of our observations. The solid grey boxes 
represent definite (SL > 4a) detections, while the black boxes indicate possible (3a < SL < 4a) 
absorbers. The W bins are 12 mA in width. 

(SL > 3a) Lya samples. Our W bins are 12 mA in width, and our upper W cutoff is 133 mA. Also 
indicated by dashed lines in Figure 6 is our least-squares fit to In [n(W)] = In [Cyv/W* ] — W/W* 
for 10 mA < W < 133 mA. For our definite sample, = 27 ± 5 mA and C w = 520 ± 179 mA. 
For our expanded sample over the same W interval, = 23 ± 3 mA and C w = 733 ± 205 mA. 

For comparison, at higher redshift, Sargent et al. (1980) obtained W* = 362 ± 21 mA, Cyy = 
154 ± 11 mA, while Young, Sargent, & Boksenberg (1982) obtained W* = 232 ± 34 mA, C w = 
135 ± 20 mA in the redshift range 1.7 < z < 3.3 for Lya absorbers with W > 160 mA. At first 
comparison, our results seem inconsistent with these results. However, as has been noted (Murdoch 
et al. 1986; Carswell et al. 1984; Atwood, Baldwin, & Carswell 1985), the parameters W* and 
Cw are highly dependent on the spectral resolution of the observations, mostly because at higher 
resolution blended lines break up into components. One sample that approximates our resolution 
and sensitivity is that of Carswell et al. (1984). Based upon 20 km s _1 resolution optical data, 
Carswell et al. (1984) obtained W* = 70 ± 20 mA for their unblended, W > 133 mA, Lya sample. 

As indicated in Murdoch et al. (1986, their Figure 3), there is evidence for a break in the slope 
of n(W) below 200 mA. As shown in Figure 7, we confirm this break, even though our statistics in 
this region (W >200 mA) are poor. For 6=25 km s" 1 , this corresponds to a break at log [A^hi]~14.0. 
The reality of this break will be discussed further in § 4, which focuses on the iVni distribution of 
the low-z Lya forest. 

At first consideration, the presence of this break is intriguing. However, one must understand 
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Fig. 6. — The differential distribution, n(W)= d 2 N 'jdz dW, plotted in natural logarithm form. 
The n(W) error bars are based upon y/~N statistics. The left panel indicates the distribution for our 
definite (SL > 4a) Lya absorbers only. The right panel includes the distribution for our expanded 
(SL > 3a) Lya sample. The W bins are 12 mA and only Lya absorbers with W < 133 mA are 
included. 



the nature of this discontinuity in relation to the Lya curve of growth to gauge its significance. As 
pointed out by Jenkins & Ostriker (1991) and Press & Rybicki (1993), this upturn in the number 
of low W clouds can be explained by the transition from saturated to unsaturated lines. The form 
of n(W)=(C\\! /W*) exp (— W/W*) is designed to model absorbers in the logarithmic or flat part of 
the curve of growth. As discussed in § 4, it has been shown that at higher column densities, Lya 
absorbers appear to follow the distribution n(Nm) = d 2 N jdz dNm = Chi-/V hi /3 . On the linear 
portion of the curve of growth, where W = ANm, 

n(W) ~ ~ d " N dNm ~ ° mN ™ or (2) 

[VV) ~ dzdw~ dz dN m dw ~ a ' or [Z) 

n(W) = W^ = A^C m W-?. (3) 

JTl. 

As derived in equation (23) for Lya, A = 54.43 x 10~ 13 mA for A^hi in 10 13 cm -2 and W in mA. 

In Figure 7 we add our absorbers with W > 133 mA to the SL > 4a In [n(W)] distribution 
of Figure 6. We increase the upper W bin size from 12 mA to 42 mA to compensate for our 
poor statistics, so that there is at least one absorber in each bin. Above W = 133 mA, we obtain 

= 259 ± 180 mA and C w = 138 ± 129 mA. Below 133 mA, we fit by equation (3), using the 
least-squares method, and find (3 = 1.89 ± 0.39 for our SL > 4a sample ((3 = 2.20 ± 0.37 for our 
SL > 3a sample) . The value of Chi is not well constrained due to the exponential nature of this 
fitting method, but the best fit gives Chi ~ 4 x 10 14 . 
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Fig. 7. — The differential distribution, n(W)= d 2 N /dz dW, of Lya absorbers with respect to W is 
plotted in natural logarithm form, for our SL > 4a Lya sample. The n(W) error bars are based 
upon y/N statistics. The W range has been expanded to include all detected Lya absorbers. Below 
W = 133 mA, the W bins are 12 mA in width; above 133 mA the bins are 42 mA wide. Horizontal 
(la) error bars are set to 1/3 of the W bin widths. 
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3.2. Integrated d/V /dz Results 

By integrating d 2 N '/dz dW from W to oo, we can determine the number density of lines per 
unit redshift, d/V /dz , stronger than W. As stated previously, we have assumed no evolution with 
z over our small range for Lya detection (0.002 < z < 0.069). Because of our very \ow-z range, 
these values for d/V /dz are a good approximation of (dAf /dz) z= o- Figure 8 shows d/V /dz , defined 

as 

POO Q2^f 

dV/d 2 (W > Wi) = ^ — d W . (4) 

The vertical axis of Figure 8 gives (d/V /dz) z= $ in terms of both W (lower axis) and log [Nm] (upper 
axis, assuming that all absorbers are single components with b- values of 25 km s _1 ). 

For comparison, Table 4 gives (dAf /dz) z= o from other HST absorption-line studies, along with 
our low-z results for their observational W limits (W m i n ). The detection limits of the various HST 
spectrographs and gratings are indicated by the dashed vertical lines in Figure 8. The unlabeled 
vertical line on the far left of Figure 8 in the GHRS/G140L+G160M merged sample of Tripp, 
Lu, & Savage (1998). In all cases, our values for (d/V /dz) z= o at various limiting W are slightly 
greater than those derived from higher limiting W studies at low-z , but are still within the la 
error ranges. The value of dAf /dz at the present epoch (z=0) is an important constraint for 
hydrodynamical cosmological models attempting to reproduce the observed baryon distribution of 
the Universe (Dave et al. 1999). The (d/V /dz) z= o values of our GHRS survey have the lowest 
measured W values to date. In § 5.1, these values will be used to place important constraints upon 
the z-evolution of dAf /dz. 



Table 4: (cbV /dz) z= o Determinations From Other Lya Studies 



Reference 


Instrument / 


z range 




(d/V/dz),= 


Our (dA//d2) z= o a 




Configuration 




(mA) 




Bit Wmin 


Bahcall et al. 1993a 


FOS b 


0< z <1.3 


320 


17.7±4.3 


18.2±6.9 


Bahcall et al. 1996 


FOS b 


0< 2 <1.3 


240 


24.3±6.6 


28.5±8.6 


Weymann et al. 1998 


FOS 


0< z <1.5 


240 


32.7±4.2 


28.5±8.6 


Impey et al. 1999 


GHRS/G140L 


0< z <0.22 


240 


38.3±5.3 


28.5±8.6 


Tripp et al. 1998 


GHRS/G140L 


0< z <0.28 


75 


71±12 


80.0±14.6 


Tripp et al. 1998 c 


GHRS/G140L 


0< z <0.28 


50 


102±16 


134.2±20.1 



"For comparison to the (djV /dz) z =o values of the previous studies, this column reports our values for (d/V /dz) z= o eval- 
uated at the minimum equivalent width limit (Wmin) of the previous studies. 
''Bahcall et al. (1993, 1996) use the G130H, G190H, and G270H HST gratings. 

c Tripp, Lu, & Savage (1998) combine their observations with observations of 3C 273 from Morris et al. (1993), who 
include observations taken with both the GHRS/G140L and GHRS/G160M. 
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Fig. 8.— Integrated <W/dz above W: dAf/de (W > W{) = f™.(d 2 N/dz dW)dW, for our 
SL > 4a sample. For clarity, we plot every third error bar. The W bins are 12 mA, connected by 
the solid line. Also shown by dotted vertical lines are the minimum W limits for the various HST 
Lya surveys presented in Table 4. The unlabeled vertical dashed line on the far left corresponds to 
the Tripp, Lu, & Savage (1998) G140L+G160M combined sample (see Table 4). The horizontal 
axes begin at W = 10 mA or log [A H i]= 12.26. 
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3.3. Rest-Frame Equivalent Width versus b Distribution 

In Figure 9, we plot the observed W distribution versus the resolution-corrected b values. In 
this figure, the 81 definite absorbers are plotted as small dots while the 30 possible (3a < SL < 4a) 
Lja detections are plotted as diamonds. In neither sample do we detect any correlation of W with 
respect to b (in the definite sample the correlation coefficient is zero). Since very broad weak features 
are indistinguishable from small continuum fluctuations, we may not be sensitive to features in the 
lower right corners (large b- values and low W) of Figure 9. A completely homogeneous distribution 
would be created in Figure 9 by adding ~15 absorbers at W < 30 mA and b = 60 — 100 km s _1 . At 
the spectral resolution of these observations these absorbers have maximum deflections of < 15% 
of the base continuum level and thus would be extremely difficult to distinguish from a slight 
undulation in the continuum. Therefore, given the absence of reliable 6-values for each absorber 
and no obvious b-W correlation, we will assume a single b- value when calculating Nm. Figure 9 
shows that there is no obvious bias in making this single 6-value choice. 



Definite Expanded 




20 40 60 80 100 20 40 60 80 100 
b (km s _1 ) b (km s" 1 ) 



Fig. 9. — Distribution of rest-frame equivalent width (W) versus 6-value of our detected Lya 
absorbers. The left panel displays the distribution of the definite (SL > 4a) sample with error 
bars. The right panel shows the SL > 4a sample as dots, combined with the 3a < SL < 4a sample 
displayed as diamonds, forming our expanded sample. The absence of absorbers in the lower right 
corners is related to an insensitivity to very broad weak features, which could be interpreted as 
continuum fluctuations. 
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4. Observed H I Column Densities 

In Table 1, we estimated the neutral hydrogen column density (Ahi) of each detected SL > 4a 
(definite) Lya absorber from its equivalent width, assuming a single component with a 6-value of 
20, 25, 30 km s -1 , and the measured 6 values. Table 2 presents iVni results for our 3a < SL < 4a 
(possible) Lya sample. As described in Paper I and in § 2, the observed 6-values have been 
corrected for the HST+GHRS/G160M instrumental profile and our pre-fit smoothing. Columns 
4 and 5 of Tables 1 and 2 indicate the measured 6- values (6 b s ) and the corrected 6- values (6), 
respectively. Figure 10 compares the results on log [./Vhi] of our sample for these various 6-values. 
Below log [A^hi]~13, all the Nm distributions of Figure 10 are similar, because the Lya absorbers 




12 13 14 15 16 
log[N HI ] (0.1 bins) 



Fig. 10. — H I column densities, ./Vhi > of our observed low-z Lya absorbers for various 6-values. 
The top three panels show the distribution based upon all features having a single 6-value of 20, 
25, or 30 km s _1 as indicated. Bottom panel shows the distribution based upon 6-values measured 
from Gaussian components and corrected by the HST+GHRS/G160M resolution element. The 
grey absorbers are our definite (SL > 4<r) sample, while the black absorbers are our possible 
(3a < SL < 4cr) sample. 
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are on the linear part of the curve of growth where Nm does not depend on the 6-value. Above 
log [A?hi]= 14, the Lya absorbers are partially saturated for 6=25± 5 km s" 1 . For Gaussian line 
profiles, the optical depth at the line center is given by: 

25kms" 1 ' 



On the logarithmic or flat-part of the curve of growth, Nm depends on the 6-values. Above 
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Fig. 11. — The \ow-z Lya number distribution per unit redshift and column density, 
n(N m )=d 2 N /dz dN m , for a constant 6-value of 25 km s _1 and log [Nm] < 14.0. The n(N m ) 
error bars are based upon y/N statistics. The left panel indicates the distribution for our defi- 
nite (SL > 4a) Lya absorbers only. The right panel includes the distribution for our expanded 
(SL > 3a) Lya sample. 

log [JVhi]~ 18.5, the Lya absorbers develop damping wings, and AT HI is once again insensitive to 
the b- value. None of our detected Lya absorbers has log [iVHi]>17. While we believe this to be a 
typical range of true values for our detections, clearly there is some uncertainty in the individual 
values of N m at log [N m ] > 13.5. 

The A^hi number density per unit redshift and column density is often modeled by a power-law 
distribution: 

In Figure 11, we display log [n(A^ni)] for both our definite (SL > 4a) and expanded (SL > 3a) Lya 
samples over the range 12.3 < log [Nm] < 14.0. Also indicated in Figure 11 are the least-squares 
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fits to log [n(Nni)] = log [Chi] — (3 log[AT HI ]. Using a 6- value of 25 km s 1 for all features, we 
obtain (3 = 1.83 ± 0.15 and log [Cm] = 12.8 ± 1.9 for our definite sample and (3 = 1.81 ± 0.12 and 
log [Chi] = 12.6 ± 1.6 for our expanded sample over the range 12.3 < log [Nm] < 14.0. There is no 
evidence for a turnover of n(Nm) below Nm = 10 125 cm -2 in either our definite or expanded Lya 
samples. The determination of (3 and Chi is insensitive to 6-value below log [Nm] < 14, since all 
absorbers are then on or near the linear portion of the Lya curve of growth, which is independent 
of 6-value. 

Above log [./Vhi]~14, the Lya absorbers become partially saturated, and the choice of 6-value 
becomes important in determining Nm- As shown in Figure 12, we detect a break in the N^f power- 
law above log [AT HI ]~14, which is possibly related to saturation or poor line statistics. However, 
Kulkarni et al. (1996) detect a similar break in the region of Nm ~10 14 ' 5 cm -2 for Lya absorbers 
at 1.7 < z < 2.1. Figure 12 compares the power-laws for the two column density regimes, 12.3 < 
log [Nm] < 14.0 and 14.0 < log [Nm] < 16.0, for both our definite (SL > 4a) and expanded samples 
(SL > 3a), again for a constant 6=25 km s _1 . Since none of the strong absorbers (log [AT HI ]> 14.0) 
are tentative detections (3a < SL < 4a), the results for the definite and expanded samples are 
identical with (3 = 1.04 ± 0.39 and log [Chi] = 1-5 ± 5.7. The large uncertainties are a reflection of 
the poor number statistics. 



To obtain better determinations of (3 and Chi , we evaluate the integrated n(Nm), 



Table 5 compares our results for both the differential, n(Nm), and integrated, I(Nm), determina- 
tions of (3 and Chi for 6-values of 20, 25, 30 km s^ 1 and 6. This table is divided into three log [Nm] 
ranges over which least-squares determinations of (3 and Chi are performed: 12.3 < log [Nm] < 14.0, 
14.0 < log [Nm] < 16-0, and the full range 12.3 < log [Nm] < 16.0. Results are given for both our 
definite and expanded Lya samples. By column Table 5 includes: (1) sub-sample name; (2) 6-value 
assumed in km s _1 ; (3-6) best-fit values for f3 and log [Chi] m cm~ 2 for both differential and inte- 
gral counts; and (7) the number of absorption lines in each subsample. Note that the results over 
the range 14.0 < log [Nm] < 16.0 are the same for the definite and expanded samples since these 
data sets are identical over this log [Nm] range. When using our measured 6-values, all features 
in this column density range occur in the same log [Nm] bin. Hence we cannot report values for 
this column density range for these 6-values. Using the integrated I(Nm) rather than the differen- 
tial distribution n(A^ni) to determine (3 and Chi is more robust since it uses the cumulative iVni 
distribution, instead of the individual Nm bin values in the parameter determination. Using the 
integrated distributions for 6=25 km s _1 , we obtain (3 = 1.72 ± 0.06 and log [Chi] = H-3 ± 0.7 for 
our definite sample and (3 = 1.81 ± 0.05 and log [Chi] = 12.4 ± 0.7 for our expanded sample over 
the range 12.3 < log [N m ] < 14.0. 
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Table 5: Differential, n(A r ni), and integrated, I(Nm), determinations of (3 and log [Chi] for the 
low-z Lya absorbers. 



12.3 < log[JV H i ] < 14.0 

Differential Integrated 



Sample 


6- value 




P 


log[C ffl j 






log[C m j 


N 


Definite 


20 


1.83 


±0.16 


12.8 


± 


2.1 


1.70 


±0.05 


11.0 


± 


0.7 


68 


Definite 


25 


1.83 


±0.15 


12.8 


± 


1.9 


1.72 


±0.06 


11.3 


± 


0.7 


71 


Definite 


30 


1.80 


±0.15 


12.5 


± 


2.0 


1.74 


±0.06 


11.5 


± 


0.7 


71 


Definite 


b 


1.84 


±0.14 


13.0 


± 


1.9 


1.77 


±0.06 


11.9 


± 


0.8 


75 


Expanded 


20 


1.81 


±0.13 


12.6 


± 


1.7 


1.78 


±0.05 


12.1 


± 


0.6 


98 


Expanded 


25 


1.81 


±0.12 


12.6 


± 


1.6 


1.81 


±0.05 


12.4 


± 


0.7 


101 


Expanded 


30 


1.80 


±0.12 


12.5 


± 


1.6 


1.83 


±0.05 


12.7 


± 


0.7 


101 


Expanded 


b 


1.83 


±0.11 


12.8 


± 


1.5 


1.86 


±0.06 


13.0 


± 


0.7 


105 



14.0 < log[AT H i ] < 16.0 

Differential Integrated 
Sample b- value /3 log [Chi J P log [Chi J ^V" 



Definite/Expanded 20 1.07 ±0.18 2.0 ± 2.7 1.34 ±0.13 6.2 ± 1.9 13" 
Definite/Expanded 25 1.04 ±0.39 1.5 ± 5.7 1.43 ±0.35 7.4 ± 5.2 10 
Definite/Expanded 30 1.40 ±0.64 6.7 ± 9.5 1.81 ±0.73 12.6 ±10.7 10 



12.3 < log[AT H i ] < 16.0 



Differential 



"log[C 



Integrated 



gc 
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Definite 

Definite 

Definite 
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Expanded 

Expanded 

Expanded 



6- value 



HI 



HI 



N 



20 
25 
30 
b 

20 
25 
30 
b 



1.44 ±0.07 
1.58 ±0.10 

1.64 ±0.12 
1.60 ±0.17 
1.50 ±0.07 

1.65 ±0.10 
1.71 ±0.12 
1.70 ±0.16 



7.6 ± 1.0 
9.6 ± 1.4 

10.3 ± 1.7 
9.8 ± 2.3 
8.5 ± 0.9 

10.4 ± 1.4 
11.2 ± 1.6 
11.1 ± 2.2 



1.53 ±0.04 
1.66 ±0.06 
1.71 ±0.08 

1.77 ±0.11 
1.57 ±0.04 
1.71 ±0.06 

1.78 ±0.08 
1.85 ±0.11 



8.9 ± 0. 
10.5 ± 
11.3 ± 1 
11.9 ± 1 
9.5 ± 0. 
11.3 ± 
12.1 ± 1 
13.0 ± 1 



81 
81 
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Fig. 12. — The differential low-z Lya number distribution per unit redshift and column density, 
n(Nm)=d 2 Af /dz dN-Ri , over the range 12.3 < log [Nyu\ < 16.0 for a constant b- value of 25 km s _1 . 
The n(Nni) error bars are based upon y/N statistics. The upper panel indicates the distribution 
for definite (SL > 4a) Lya absorbers only. The lower panel includes the distribution for expanded 
(SL > 3a) Lya sample. Indicated on each panel are the best-fit power-laws for each log [Nm] 
regime, 12.3 < log [Nm] < 14.0 and 14.0 < log [iVm] < 16.0. In the lower JVhi regime, the log [Nm] 
bins are 0.12, while in the upper regime the bins are 0.40. 
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As indicated in Table 5, the values for (3 that we obtain over the column density range 12.3 < 
log [Nm] < 14.0 are consistently >1.7. This value is in agreement with our results of § 3, where we 
fitted n(W)(x W' 13 and obtained (3 = 1.89 ± 0.39. We derive our best values from the integrated 
SL > 4a sample, using a constant value 6=25 km s _1 : 

P = 1.72 ± 0.06 , log [C m ] = 11.3 ± 0.7 for 12.5 < log [N m ] < 14.0, (8) 

(3 = 1.43 ± 0.35 , log [C ffl ] = 7.4 ± 5.2 for 14.0 < log [N m ] < 16.0 , (9) 

(3 = 1.66 ± 0.06 , log [Cm] = 10.5 ± 0.8 for 12.5 < log [N m ] < 16.0. (10) 

These results are in contrast with the recent higher-z results of Kim et al. (1997), Lu et al. (1996), 
and Hu et al. (1995), who find that /3~1.4 adequately describes n(Nm) for 12.3<log [./Vhi]<14.3 over 
the redshift range 2.17<z<4.00. In conjunction with other studies between z=4 and z=2, these 
authors claim that the n(Nui) break appears to strengthen and move to lower column densities with 
decreasing redshift. Furthermore, Kim et al. (1997) report a break in n(Nm) above log [Njji] > 
14.3 that steepens {(3 = 1.7 — 1.8) and moves to lower column density with decreasing redshift 
(2.17<z<4.0). In contrast, at low-z we see a flattening in n(Nm) above log [A^hi]=14.0. It is 
possible that the apparent flattening in our data is an artifact of our selection of a constant b- 
value for all Lya absorbers. If the stronger Lya absorbers had larger 6-values, then the inferred 
column densities of these features could be much smaller. This would leave us with fewer features 
above log [AT HI ]>14, and therefore no evidence for a break in n(Njn) at log [Nm] = 14.0. In § 5.3, 
we will examine the possibility that the break in n(Nm) is a consequence of the z evolution of 
d 2 N/dz dNm ■ 
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5. Observed Redshift Distribution 

In this section, we examine the redshift distribution of the low-z Lya forest. In particular, we 
want to examine the evidence for structure at specific redshifts, or for any evolution of d/V /dz with 
z. We begin with Figure 13, a presentation of the observed number distribution per redshift bin 
without any sensitivity correction, M(z), for both our definite and expanded Lya samples. The 
left axis corresponds to the displayed histograms oiJ\f(z). Our redshift range, 0.002 < z < 0.069, 
is divided into 12 redshift bins, 5z= 0.0056. Individual SL > 4a Lya absorbers are plotted as 
pluses, whose W(mA) is shown on the right vertical axis. Individual 3a < SL < 4a Lya absorbers 
are plotted as 'x's. Like the Af(W) distribution, the N(z) distribution is not the true redshift 
distribution, because we have not corrected for the varying wavelength and sensitivity coverage of 
our observations. 

We show the correction for our varying wavelength and sensitivity coverage, Az(z), in Fig- 
ure 14. In this figure, the panels show various stages of correction for the cumulative distribution 
of Az(z) for SL > 4a detections in all sightlines. The upper panel of Figure 14 presents the full 
pathlength availability of our sample. The left axis gives cAz(z) in Mm s _1 , while the right axis is 
in redshift units. Subsequent panels present the available pathlength after removing portions of our 
wavelength coverage due to obscuration by extragalactic non-Lya features and intrinsic features 
(Intrinsic+Non-Lya), by Galactic+HVC features, and by our cz cm - 1,200 km s" 1 proximity limit. 
The bottom panel presents the pathlength available after removing all spectral regions not suit- 
able for detecting intervening Lya absorbers. The available pathlength in Figure 14 approximately 
represents Az(z) for features with W>150 mA. As presented in Figure 4, our pathlength is ap- 
proximately constant at cAz(W)=0.387 for W>150 mA. To properly characterize the low-z Lya 
absorber distribution in redshift, one must accurately account for varying available pathlength as a 
function of both 7V HI and z of the observations. Figure 15 displays the combined two-dimensional 
sensitivity function, cAz(Nm,z), for GHRS/160M observations as a function of z and Ahi, after 
accounting for the aforementioned corrections and obscurations. 

Unfortunately, we do not have enough statistics (absorbers) to fully analyze the Nm versus 
z distribution over the small z-range of our sample. However, there is no obvious trend over our 
local pathlength. We therefore assume no z evolution (Weymann et al. 1998) of the column density 
distribution for our Lya sample {(3 is constant for all z in our observed range). 

It is customary to describe the pathlength-corrected number dependence of Lya absorbers as 
a power-law in redshift, n(z) = C z (l + z) 1 . Because our spectra vary in wavelength coverage and 
Ahi sensitivity, we must correct n(z) for incompleteness. In particular, we know that there exists 
a distribution, 

^ = A ?N N ?AN = °™ ■ (ID 
Az(N m )AN m 

that describes the Ahi distribution, and we know the pathlength availability, Az(Nm,z), as a 
function of both z and Nm. Therefore, the z distribution evaluated at each redshift bin, n(zi), can 
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Fig. 13. — Histogram of A^(^) and distribution of rest-frame equivalent width (W) versus redshift. 
The histograms correspond to the left axis and indicate M(z), the number of absorbers per redshift 
bin (Sz). The solid histogram indicates the distribution of the definite (SL > 4a) Lya sample, 
while the dotted distribution is for the expanded (SL > 3a) Lya sample. The pluses ('+') and x's 
correspond to the right axis and indicate W of the individual Lya absorbers in the definite and 
possible (3a < SL < 4a) samples, respectively. 
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Fig. 14. — Cumulative available pathlength (Az) for all sightlines as a function of z, for features 
with W > 150 mA. Left axis gives cAz(z) in units of Mm s _1 , while the right axis gives Az(z) in 
terms of z. The upper panel presents the full, uncorrected Az(z). Subsequent panels present Az(z) 
after removing portions of our wavelength coverage due to obscuration by extragalactic non-Lya 
features, Galactic + HVC features, and the (cz em - 1,200 km s _1 ) proximity limit. The bottom 
panel presents our cAz(z) for detecting intervening Lya absorbers after applying these corrections. 
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Fig. 15. — Bivariate pathlength sensitivity correction, cAz(Nui,z). The sensitivity function corrects 
d 2 M/dz dNni by accounting for varying regions of Nm sensitivity and z coverage. The vertical 
axis indicates the available pathlength (in Mm s" 1 ) for a given column density, Nm, and redshift 
(wavelength). When viewed from the column density axis, this figure reduces to the dotted line in 
Figure 4, and when viewed from the redshift axis, the bottom plot of Figure 14 is recovered. 
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be corrected for incompleteness by, 

n(N m )n( Zi )dN m = / J - ^ HI ' lj diY ffl , (12) 

in ■ Jn ■ ozofsm 

Iv miii JY min 



f 

Jn, 



f 

JN, 



-dN m (13) 



l Ninin Az(N m , Zi )AN m (N m ) 
or, 

r°° N"(Nm,zi) j at 

( N _ Mnin A Z (Af H I, Zi )AA r Hl(A f Hl) aiVHI /-.^ 

ijv min n ( N ™)dN m 

where, 

n(AT m ) = Cm iV H f- ( 15 ) 

The separability of the bivariate distribution, d 2 N/dz dNm, assumed in equations (12) and (14) 
implies that there is no z evolution of the Nm distribution of our Lya absorbers within the small Az 
of our GHRS observations. Under this assumption, the integrals in equation (14) can be combined, 
and the observed redshift distribution of absorbers, n(z), can be expressed as: 

n(zi)= I" jV(N m ,Zi) =CJ1 + Zi y. (16) 

JN mill CmN^Az(N m ,z i )AN m (N m ) 

Our goal is to search for the variations in dj\f /dz with redshift, accounting for the sensitivity limits 
of our spectra. Thus, n(zi) in Figure 16 shows dj\f /dz(zi) normalized by the expected oW /dz given 
our observed dAf / dNm distribution and sensitivity limits for a uniform distribution of clouds over 
our redshift range. For all z bins, -/V m i n is set to 10 12 ' 3 cm" 2 , although most bins have no pathlength 
at this column density. 

The top panel of Figure 16 displays n{z\) for the SL > Act sample. Error bars in this figure 
are based upon Poisson statistics. The increase in features at z=0. 055-0. 065 is due to a single 
cluster of lines in the upper portion of the PKS 2155-304 spectrum. We specifically selected this 
sightline to observe this complex of lines, therefore, we trimmed this portion of the PKS 2155- 
304 sightline to remove any bias introduced by selecting a sightline which we previously knew 
contained an strong complex of Lya absorbers. The trimmed distribution is shown in the lower 
panel of Figure 16. Above z=0.05, we have very little pathlength after removing the upper portion 
of the PKS 2155-304 spectrum from our sample. The large uncertainty in n{zi) at z=0.04 is due 
our reduced pathlength in this bin due to the presence of Galactic S II A1259 and Si II A1260.4 + 
Fe II A 1260. 5. We see no compelling evidence for any z evolution over our small range in redshift 
(7 is consistent with 0; 7 = -0.01 ± 0.02). 



5.1. Redshift Evolution of dM/dz. 



Figure 17 displays dN /dz over the redshift interval < z < 3, for several studies over two 
W ranges: W>240 ml (log [N m ]>U for 6=25 km s" 1 ) and 60 mA < W < 240 mA (13.1 < 
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Fig. 16. — The distribution, n(zi) or normalized cUV /dz per unit redshift, is plotted versus z for the 
definite (SL > 4a) sample in the upper panel, and for the trimmed (SL > 4a) sample in the lower 
panel. Error bars are based upon Poisson statistics. The excess in n(zi) above z=0.05 in the upper 
panel is due to a complex of lines in the PKS 2155-304 sightline at A > 1280 A. We specifically 
selected this sightline to observe this complex of lines, therefore, we trimmed this portion of the 
PKS 2155-304 sightline to remove any bias introduced by selecting a sightline, which we previously 
knew contained an strong complex of Lya absorbers. No obvious z evolution of dM /dz in observed 
over our small redshift range. The bin size is (5z=0.0056. For comparison, a value of n(z)=0.5 
corresponds to dAf /dz ~90, integrated over all Nm to which we are sensitive. 
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log [Njjx] < 14.0 for 6=25 km s r ). The lower distribution in Figure 17 is data normalized to 
W > 240 ml by Weymann et al. (1998), while the upper distribution corresponds to absorbers 
in the range 13.1 < log [Nm] < 14.0. The data points indicated by squares (log [1 + z]<0.4) were 
obtained as part of the HST/FOS Key Project (Weymann et al. 1998), and while normalized to 
W > 240 mA, does contain a few absorbers below that limit, down to 60 mA. The stars and 
diamonds correspond to ground-based data (log [1 + z]>0.4), taken with an equivalent width limit 
of W > 360 mA, reported by Lu, Wolfe, & Turnshek (1991) and Bechtold et al. (1994). Since 
the distribution of oW /dz at W > 240 mA can be described by a single power-law in z, the high- 
z data were scaled by Weymann et al. (1998) to be consistent with the HST/FOS Key Project 
W > 240 mA data for z < 1.5. The two low- z points indicated by filled circles are taken from this 
survey, one point for each of the two sensitivity ranges. The solid lines indicate the best fits to 
the W>240 mA data above and below z= 0.4, and have slopes of 7 = 0.16 ± 0.16 (log [1 + z]<0.4, 
Weymann et al. 1998) and 7 = 1.85±0.27 (log [1 + z]>QA, Bechtold et al. 1994). The solid triangle 
represents the mean of the Kim et al. (1997) data (open triangles), for which 7 = 1.19±0.45 (dash- 
dot line). While the difference between 7 = 1.85 for W > 240 mA and 7 = 1.19 for lower- W is 
suggestive of a slower evolution at high-z for the lower W absorbers, the large error bars on the 
Kim et al. (1997) data mean that this is far from definitive at this time. However, the Savaglio 
et al. (1999) point, using the QSO in the Hubble Deep Field-South, is suggestive of a somewhat 
slower evolution as well, but again with large error bars. 

The interpretation of the high-W data follows the numerical modeling of Dave et al. (1999). 
The break in dN /dz at z~0.4 can be explained by the expansion of the Universe, combined 
with a rapid decline in the metagalactic ionizing flux at z<2. The expansion of the Universe 
causes Lya absorbers to decrease in density with decreasing redshift in the absence of gravitational 
confinement, resulting in a rapid decrease in recombinations and thus in the observed column 
density of H I for any initial baryon overdensity. Also, since there exists an inverse relationship 
between the number of clouds and Nm, as the Universe expands there are fewer clouds at any 
given W. At z<2, the rapid decrease (Shull et al. 1999b) in the ionizing background intensity, 
Jv(z), allows Lya clouds to become less ionized (whi oc n^/Jo), which results in an increase in the 
amount of hydrogen detectable in Lya. In other words, at z<2, the decline in the recombination 
rate caused by the expansion of the Universe is countered by a decreasing photoionization rate (Fhi) 
due to the declining UV background, resulting in the dramatic break in dAf /dz seen in Figure 17 
at log[l + z\ ~ 0.4. 



5.2. An Important Technical Issue for Redshift Evolution of low W Absorbers 

It has been suggested by cosmological hydrodynamical simulations (Dave et al. 1999) that one 
would expect lower-W Lya absorbers to evolve more slowly (smaller 7). Evidence for this effect 
is seen in Figure 17 for log [1 + z]>0.4 in that the lower-W distribution has a shallower slope (7= 
1.19; Kim et al. 1997) than the higher-W distribution (7=1.85; Bechtold et al. 1994). While 
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this difference is suggestive of a slower evolution at high-z for the lower W absorbers, the large 
error bars on the Kim et al. (1997) and Savaglio et al. (1999) data mean that this result is not yet 
definitive. A slower evolution for lower W absorbers has also been reported in HST/FOS data at 
lower redshift by Weymann et al. (1998). Indeed, the best-fit value for the lowest- W absorbers in 
the Key Project data has a very slight negative evolution (i.e., fewer absorbers at higher-z). To 
test this hypothesis for the low-W absorbers with 13.1 < log [iV HI ] < 14.0, we connect the mean 
point of Kim et al. (1997, upper dashed line) to our z=0 data point (large filled circle in Figure 8) 
giving a slope 7 = 0.75 ± 0.15. A similiar analysis for the W > 240 mA sample (connecting the 
mean point of Bechtold et al. (1994) to the extrapolation of the Weymann et al. (1998) best-fit line 
to z=0) produces a nearly identical result, 7 = 0.75 ± 0.09, indicating no distinguishable difference 
in the overall z-evolution of cbV /dz for these two distributions. In some ways, the much larger 
redshift difference between these data points yields a much more secure result than either the Kim 
et al. (1997) or Weymann et al. (1998) analyses alone and also one that apparently differs from 
these earlier papers. But because these other works are measuring only portions of the redshift 
evolution, while our two-point measurement is for the full z range, these results can differ and still 
not be contradictory (i.e., the "breakpoint" from faster to slower evolution could occur at lower 
redshifts for the low-W absorbers). 

Nevertheless, we suspect that the source of this apparent discrepancy may arise from the 
manner in which these earlier analyses have fitted the observed line densities. Specifically, Dave 
et al. (1999) have extended the n(W) = (Cw/W*) exp (— W/W* ) relationship to inappropriately 
low W values. As discussed in § 3.1 and as shown in Figure 7, below Ws of about 133 mA, where 
the absorption lines are on the linear portion of the curve of growth, the n(W) relationship is 
better expressed as n(W) oc W _/3 . Thus, the Dave et al. (1999) fitting technique results in an 
overestimate of W* and an underestimate of 7 by assuming an exponential distribution of absorber 
numbers rather than the steep power-law, which is observed. While this overestimate almost 
certainly occurs in the Dave et al. (1999) analysis, the HST/FOS Key Project line list does not 
contain enough absorbers at W < 150 mA to cause a significant overestimate by assuming an 
exponential distribution in W. The fact that Weymann et al. (1998) find similiar W* values for all 
their analyzed subsamples is suggestive that this is, in fact, the case. On the other hand, the Dave 
et al. (1999) synthetic line list almost certainly has many absorbers that are on the linear portion 
of the curve-of-growth. This would explain why the Dave et al. (1999) analysis of their simulations 
significantly overpredicts (factor of 2.3) the number of W > 50 mA Lya absorbers found by Tripp, 
Lu, h Savage (1998), and would presumably overpredict the number density of low W absorbers 
that we have found as well. 

The amount of this underestimate of 7 increases with decreasing redshift in the Dave et al. 
(1999) analysis (their Figure 8), casting doubt on their prediction that weaker W features should 
evolve faster in dN /dz. Because this overestimate of the line density and attendant underestimate 
of 7 are most important at W < 200 mA, it is the lowest W point in Figure 9 of Dave et al. (1999), 
and possibly in Figure 7 of Weymann et al. (1998) that are too low. At higher W, where the inferred 
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equivalent widths are not affected by this technical issue, both the Dave et al. (1999) simulated 
data and Weymann et al. (1998) Key Project data still show marginal evidence for slower evolution. 
Thus, if the proper analysis were made of both the simulations and the local Lya data, they would 
still agree. However, there would be little evidence that the higher and lower W absorbers evolve 
at different rates between z=3 and 0. This is exactly what our data show in Figure 17. When our 
data are combined with that of Kim et al. (1997) and Weymann et al. (1998), there is no conclusive 
evidence for a difference in evolutionary rates between high- and low-W absorbers. 

The lower W results are complicated by the lack of dAf /dz data in the range 13.1 < log [Nm] < 
14.0 for 0.1 < log [1 + z] < 0.4. This redshift range (0.25 < z < 1.5) is indicated in Figure 17 by 
the question marks. Two HST+STIS cycles 8 & 9 projects (B. Jannuzi, PI) are scheduled to obtain 
spectra of sightlines in the redshift range, 0.9 < z < 1.5 to address this deficiency. The improved 
sensitivity of HST+COS also should be able to provide data to clarify the evolution of weak Lya 
lines at redshifts < z < 1.5. Furthermore, the inclusion of our cycle 7 STIS observations (J. T. 
Stocke, PI) should double our number statistics at low redshift and help determine the low-W line 
density more precisely as well. 

5.3. Redshift Evolution of d 2 M/dz dN m . 

In this section, we examine the redshift evolution of the bivariate number distribution with 
respect to redshift and column density, d 2 N ' /dz dNm- If any redshift evolution is detected, it 
could yield insight into the merging or dissipation of low-z Lya clouds. In Figure 18, we plot 
d 2 M /dz dNm, multiplied by 7V"hi to expand the structure near log [A^hi] ~ 14. We only include 
our results below log [Nm] < 14, since above this column density value our statistics become poor. 
However, we combine our log [Nm] < 14 results (z~0) with the HST+FOS Key project data of 
Weymann et al. (1998) in the redshift range < z < 1.3. Note that there is good agreement 
between the Weymann et al. (1998) lowest iVm points and the highest A^hi points from our survey. 
For comparison, Figure 18 also includes data at higher redshift, (z) = 3, compiled by Fardal, 
Giroux, & Shull (1998). 

The preliminary suggestion from Figure 18 is that the distribution of Lya absorbers moves 
leftward and downward from z=3 to the present. The major factors governing the z evolution of 
dAf /dz for the Lya forest are: (1) the expansion of the Universe; (2) the rapid decline in J u {z) at 
z<2; and (3) the merger and dissipation of Lya absorbers. The expansion of the Universe will tend 
to disperse clouds, causing d 2 N ' /dz dNm to move leftward in Figure 18. In addition, the expansion 
is expected to reduce d 2 M /dz dNm a t any given A^hi, since fewer clouds are available at higher 
column densities, which have been reduced in column density down to A^hi- These are the same 
effects that drive dM /dz to lower values above z>2 (Figure 17). At z<2, the rapid decline in the 
ionizing background, J u (z), will cause the ionized fraction in the Lya absorbers to drop. This will 
increase the measured Nm, countering the effect of expansion and causing dM /dNm to move to 
the right in Figure 18. These same two factors cause the break in cbV /dz (Figure 17) for the higher 
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Fig. 17. — Comparison of log [dM /dz] versus log [1 + z] for two iVni ranges. The lower distribution 
corresponds to W > 240 mA (A%i > 10 14 cm -2 for 5=25 km s _1 ). The upper distribution 
corresponds to absorbers in the range 13.1 < log [AT HI ] < 14.0. The z~0 points (solid circles) 
are taken from our survey (Figure 8), for each of the two iVni ranges. Solid lines are taken from 
Weymann et al. (1998) and have slopes of 7 = 0.16 (log [1 + z]<0A) and 1.85 (log [1 + z]>0A). 
The filled triangle is the mean of the Kim et al. (1997) data (open triangles), for which 7 = 1.19. 
Connecting this point (upper dotted line) to our datapoint with 13.1 < log [A^hi] < 14.0 at z~0 
gives 7 = 0.75 ± 0.15. For the W > 240 mA sample of Weymann et al. (1998), the nearly identical 
result of 7 = 0.75 ± 0.09 is obtained, indicating no difference in the overall z-evolution of cW /dz 
for these distributions. The complete evolutionary picture for low-A%i absorbers is not available 
due to the lack of data for 0.1 < log [1 + z] < 0.4, indicated by the question marks. 
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column density absorbers. These two factors alone can explain the offset in Figure 18 between 
the lower and higher redshift d 2 N /dz dNm distributions. Additionally, we expect that some Lya 
absorbers are merging into higher column density systems, or collapsing due to gravity. This will 
cause additional changes in the Nm structure of d 2 J\l ' /dz dNm- This evolutionary track should be 
considered preliminary, as more comprehensive data sets, with better statistics over wider column 
density and redshift ranges, are needed to explore the evolution of d 2 N jdz dNm in detail. 



5.4. The Opacity of the Low-z Lya Forest 

For Poisson-distributed Lya clouds (Paresce, McKee, & Bowyer 1980), the effective continuum 
opacity is given by 

T c s{v bs,z obs ,z) = J dz gjv ffl q z t 1 ~ ex P(- T H)] dN ni , (17) 

where d 2 N /dNm dz is the bivariate distribution of Lya absorbers in column density and redshift, 
and t{v) = Ahi chi(^) is the photoelectric (Lyman continuum) optical depth at frequency v = 
^o6s/(l + z ) through an individual absorber with column density iVm- For purposes of assessing 
the local attenuation length, it is useful (Fardal & Shull 1993) to use the differential form of 
equation (17), marking the rate of change of optical depth with redshift, 

dT cS f°° d 2 M 



dz 



r°° d z N 

L TmnrJ 1 - exp{ - T{,y))]dNm ■ <18> 



The attenuation length, in redshift units, is then given by the reciprocal of dr e s/dz. For simplicity, 
we calculate dr e g/dz at the Lyman edge (912 A). This is a reasonable approximation due to the 
strong dependence of c"hi(A) on A. As dr e g/dz approaches and surpasses 1, it significantly affects 
the radiative transfer of the metagalactic ionizing background. 

As shown in Figure 19, the cumulative opacity of low- z Lya clouds is dr e ff/dz~0.01 for 
log[iVHi]< 13, rising to ~0.1 for log[A^ni]< 15. Figure 19 gives dT c s/dz for three 6-values, 20, 25 
and 30 km s _1 for all Lya absorbers. The crossover with respect to 6- value in the dr e g/dz curves 
between 14.5 < log[iVHi]< 15.5 arises from the Lya curve of growth and small number statis- 
tics. A lower assumed value of b will systematically increase the inferred column density for lines 
above log [A^hi]~14. For our limited data set, this reduces the observed number of absorbers with 
log [./Vhi]=14.5-15.5 for 6=20 and 25 km s -1 . Compared to 6=30 km s _1 , this causes a deficiency 
in dr e ff/dz for lower 6-values in this limited column density range. Once all absorbers are ac- 
counted for, the cumulative dr e g/dz becomes larger for lower assumed 6-values. At log [Ahi]>15, 
dTeff/dz becomes uncertain due to both poor number statistics in our sample and to the large 
dependence on 6-value. For a constant 6=25 km s _1 for all low-z Lya absorbers, dr e ff/d2;~ 0.2 for 
log [A^hi]< 16. However, if 6=20 km s _1 is a more representative Doppler parameter, dr e fj/dz~0.4 
for log [Ahi]< 16. If, as inferred from higher redshift studies (Hu et al. 1995) and from ORFEUS 
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Fig. 18. — Comparison of the column-density distribution, Nm x (d 2 J\f/dz dNnij, for z=0 (our 
data), for the HST/FOS Key Project at 0<z<1.3 (Weymann et al. 1998) and for (z}=3. For our 
data at (z)=0, we used a constant value of 6=30 km s _1 in determining the column densities to be 
consistent with the FOS data. We include only our results below log [iVHi]<14, since our statistics 
become poor at higher Nw The (z)=3 data are taken from Fardal, Giroux, & Shull (1998) and 
references therein. The solid lines connect the displayed data with a two-point summation, and are 
indicative only. 
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and FUSE Ly/3 data at low- z (Hurwitz et al. 1998; Shull et al. 2000), some Lya clouds have 
6~15 km s _1 , the cumulative Lya cloud opacity in the local Universe could approach unity for 
log [A^hi]~16. As Figure 19 indicates, Lya absorbers with 15<log [7vhi]<17 probably dominate 
the continuum opacity of the low- z Lya forest and could impact the level of the ionizing back- 
ground (Shull et al. 1999b). Characterizing the distribution of these absorbers accurately at low- z 
will remain a challenge, even for HST+COS, but would be very important in understanding the 
extragalactic ionizing background in the current epoch. 




Fig. 19. — Cumulative opacity, d.T e g/dz, at low redshift at the Lyman edge (912 A) as a function 
of Nm, calculated for three different 6- values: 20, 25, and 30 km s _1 . The crossover of the curves 
for dr e ff/dz between 14.5<log [A^hi]<15.5 arises from the Lya curve of growth and small number 
statistics (see text). 
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6. Lya Absorber Linear Two-Point Correlation Function 

The two-point correlation function (TPCF, £) can be estimated from the pair counts of Lya 
absorption lines along each line of sight in our data according to : 

Here, N Q ^ S is the number of observed pairs and iV ran is the number of pairs that would be ex- 
pected randomly in the absence of clustering, in a given velocity difference bin, Av. We determine 
iV ran from Monte Carlo simulations based upon our determined number density, d 2 N /dz dNm , 
as well as the wavelength extent and sensitivity limit of our observations. Like the pathlength 
normalization vector, we include only those portions of the spectra not obscured by Galactic lines, 
non-Lya lines, and spectral regions blueward of cz ern - 1,200 km s _1 of the target. 



At each position along the spectrum, the probability of finding an absorber is calculated by: 

'N min (\) oz dN m JN min {\) 



P ( A ) = / / x T^Tt dN ™ dz ~ Az W / Chi N~f dN m , (20) 



where N m i n (X) is based upon the sensitivity limit of the spectrum. The integral in z can be 
replaced by the z width of each pixel, Az(A), since there appears to be no z evolution between 
0.002 < z < 0.069 (i.e., 7 = 0). The quantities (3 and Chi were taken from our expanded sample 
over 12.3 < log [N m ] < 16.0 for 6=25 km s~\ and have values of 1.71 ± 0.06 and 11.3 ± 0.6, 
respectively. The probability, -P(A), is then compared to a uniformly distributed random number. If 
the probability exceeds the random number, an absorber is inserted into the Monte Carlo simulation 
at this position (A). To correct for blending effects, once an absorber is inserted into the Monte 
Carlo simulation, P(X) is set to zero for the adjacent 12 pixels. This corresponds to 2.5 resolution 
elements or ~50 km s _1 , since no pairs were observed at our resolution with separations less than 
50 km s _1 . Undoubtedly, such closer pairs exist, but at our resolution we are insensitive to them. 
Because P(\) depends exclusively on our observed distribution, any blended lines counted as a 
single absorption system will affect AT obs and iV ran identically, leaving £(Av) unchanged. The mean 
b- value for our combined pre- and post-COSTAR SL > 4a sample is ~38 km s -1 , corresponding 
to a Gaussian width (Wg) of 27 km s _1 . One would expect to begin resolving pairs separated 
by two Gaussian widths, which is in agreement with our observed 50 km s _1 cutoff. Our pre- 
COSTAR b- values have a higher median value of ~60 km s _1 , or a Gaussian width Wg~42 km s _1 . 
Therefore, in our Monte Carlo simulations, we may slightly overestimate the number of random 
pairs < 70 km s~ x in the pre-COSTAR sample, or underestimate £ in our lowest velocity bin. For 
each sightline, we performed 1,000 simulations (N s ) and combined them to form A" ra n- The error 
in N ran (Av), denoted a ran , is taken to be y 7 'N T!m (Av) . For proper scaling of £(Au), both iV ran (Af) 
and a ran are normalized by N s . 

Table 6 lists all absorber pairs with velocity separations of 50 < Av < 150 km s _1 , and 
Figures 20 and 21 display the results of our TPCF analysis, £(Au). Table 6 lists by column: (1) 



-42 - 



The central wavelength of the line pair; (2-3) the wavelength and rest-frame velocity separation of 
the pair; (4-5) the equivalent widths of the two absorbers; (6-7) the observed 6-values of the two 
absorbers; and (8) the target sightline. For this analysis we have excluded one Lya feature (and 
0.4 A of pathlength) at 1226.96 A in the PKS 2155-304 sightline because it could be Si III A1206.5 
absorption related to the strong Lya line at 1236.43 A. This line would have combined with 
another weak Lya line to produce a pair in Table 6. We have visually inspected the other line 
pairs in Table 6 and find that all but two entries clearly are pairs of distinct lines. The remaining 
two pairs (Fairall 9 and the A cen =1236.21 A pair in PKS 2155-304) could be broader lines which 
have been subdivided by our profile fitting routines. Higher resolution spectra are required to be 
certain. But, at most, the TPCF peak at Av < 150 km s^ 1 may be overestimated by 15% (2/13). 
We have also removed the A > 1280 A portion of the PKS 2155-304 sightline owing to a strong 
cluster of lines (Shull et al. 1998). Based upon lower resolution data (Bruhweiler et al. 1993), the 
PKS 2155-304 spectrum was obtained specifically to study this cluster of absorbers around 1290 A. 
Because of these special circumstances, we also exclude this portion of our sample from our TPCF 
analysis. 

The distribution drawn in Figure 20 with a solid line in the upper panel displays the number 
of observed Lya pairs with the indicated velocity separations (Av ) uncorrected for the varying 
wavelength coverage and spectral availability of our observations. The velocity separations between 
any two absorbers along the same line of sight are calculated by: 

A -T^r <*> 

where Az = Z2 — z\ = (A2 — Ai)/1215.67 A and (z)= (z2 + z\)/2. One concern is that we might have 
misinterpreted weak metal lines as Lya absorbers, as mentioned above. The vertical dotted lines in 
Figure 20 indicate the rest-frame separations between Lya and expected metal lines (Si III A1206.5, 
N V AA1238, 1242, S II A1250, or S II A1253), indicating that this is not a concern. The dotted line 
is the random distribution, N ran (Av), which accounts for the varying sensitivity and wavelength 
coverage of our observations and leads to £(Au), which is displayed in the bottom panel. The 
velocity separation bins in Figure 20 are cAz= 100 km s -1 . 

Figure 20 shows an excess in £ at the lowest velocities, indicating clustering in the local Lya 
forest, and broad deficits centered at 1850 and 3500 km s _1 . Other features at larger separations 
are also observed, but are less than 2a. Figure 21 was constructed with bin size Av= 70 km s^ 1 to 
show the effects of different binning of £(z). Calculating the significance of any peak or deficit in 
the TPCF is achieved by summing £/<7£. We find a 3.6cr excess in the velocity range 50-250 km s _1 
in the 100 km s _1 bins, £(50 — 150 km s _1 ) = 1.9 ± 0.5, and a 4.5cr excess in the velocity range 
50-260 km s _1 with the 70 km s~ x bins. 

For comparison, Ulmer (1996) calculated the TPCF for Lya lines with W > 240 mA in the 
range < z < 1.3 using data from the HST/FOS Key Project (Bahcall et al. 1993, 1996). He 
concluded that the W > 240 mA local Lya forest is clustered on the scale of Av < 500 km s -1 , 
with amplitude £(250-500 km s" 1 ) = l.&t{ m %- The deficit in £(Av) at Av~l,700 km s" 1 is also 
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Fig. 20. — Two-point correlation function (TPCF, £) of the Lya absorbers as a function of velocity 
pair separation, Av, in 100 km s _1 bins. Upper panel (solid line) gives the number of observed 
Lya pairs, N ^ s (Av), uncorrected for the varying wavelength coverage and spectral availability of 
our observations. The vertical dot-dashed lines are the Av separations between Lya and nearby 
metal lines (Si III A1206.5 and N V AA1238, 1242 are indicated, the unlabeled lines are S II). 
The dotted line is the random distribution, iV ran (Ai;), accounting for the varying available spectral 
regions. Bottom panel gives the corrected TPCF, £(Ai>). The 1 and 2a ranges of are indicated 
by the short and long dashed lines, respectively, in the bottom panel. The A > 1280 A and 
1226.76 < A < 1227.16 A regions of the PKS 2155-304 spectrum have been removed. 
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Fig. 21. — Two-point correlation function as explained in Figure 20, with velocity bins of 70 km s 1 ; 
each bin is 1 Mpc. 
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present in the FOS data, but is unexplained by Ulmer (1996). Interestingly, when metal-line 
systems are included in Ulmer's analysis, this deficit mostly vanishes. At high redshifts, there have 
been numerous reports of marginal detections of excess power in the TPCF at Av < 150 km s _1 
(Chernomordik 1995; Hu et al. 1995; Kulkarni et al. 1996; Lu et al. 1996), as well as reports 
of non-detections over the same scales (Pettini et al. 1990; Rauch et al. 1992). Impey, Petry, & 
Flint (1999) recently constructed a TPCF for z~0 Lya lines that showed no excess power at 
Av > 500 km s" 1 using GHRS/G140L data. Their TPCF was not sensitive to smaller velocity 
separations, owing to the lower resolution of the G140L. The Act peaks in our low-z TPCF figures 
at Af~100 km s~ x are at the same Av as the 2-3<r peaks at high redshift. Because this peak 
occurs at small separations, close to our resolution limit, it is the least likely feature to be caused 
by the limited wavelength coverage of our spectra. It is also unrelated to the velocity bin size, 
since different bin sizes do not alter its significance appreciably. Since no other low-z study has 
the velocity resolution required to see this peak clearly, the line pairs in Table 6 would not be fully 
resolved by HST/GHRS G140L or FOS spectra like those used by Impey, Petry, & Flint (1999) 
and Weymann et al. (1998). Therefore, we believe that the TPCF of the Lya forest at low-z has 
the same general characteristics as at high-z, but with a slightly greater excess at low Av. 

Due to the rather complex detection limit, and low pathlength, of our spectra, we are less 
confident of the reality of the possible excesses and deficits above Av=500 km s _1 . These features 
do not appear to be present in either high-z (Rauch et al. 1992) or low-z (Impey, Petry, & Flint 
1999) studies. However, if these features are real, they may be related to the local widths of 
filaments (Av=500 km s _1 ) and galaxy voids (Av=600 - 1,900 km s _1 ), since there is a possible 
excess of Lya clouds associated with filaments compared to voids (Stocke et al. 1995; Shull, Stocke, 
& Penton 1996). 

In HST cycle 7, our group has observed 13 sightlines with the STIS/G140M, which has reso- 
lution comparable to GHRS/G160M. These data should approximately double the number of Lya 
absorbers in our sample, and thus will provide much improved information on the low-z lya forest 
TPCF. We will report on the full-sample TPCF in a later paper (Penton, Stocke, & Shull 2000c, 
Paper IV). 
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Table 6: Lya Line Pairs with Velocity Separations of 50-150 km s 





AA 


AV 


Wi 


W 2 


h 


b 2 


Sightline 


(A) 


(mA) 


(km s" 1 ) 


(mA) 


(mA) 


(km s _1 ) 


(km s _1 ) 




1247.76 


0.38 


88 


20 


146 


37 


31 


AKN120 


1247.88 


0.62 


145 


20 


64 


37 


28 


AKN120 


1248.07 


0.24 


57 


146 


64 


31 


28 


AKN120 


1264.89 


0.42 


95 


30 


27 


46 


28 


FAIRALL9 


1247.76 


0.35 


82 


39 


67 


25 


36 


H1821+643 


1247.37 


0.32 


74 


48 


21 


32 


37 


MARK279 


1241.66 


0.30 


69 


57 


39 


24 


26 


MARK279 


1236.60 


0.60 


142 


207 


24 


85 


28 


MARK817 


1235.87 


0.25 


59 


64 


81 


69 


60 


PKS2155-304 


1236.21 


0.43 


102 


81 


218 


60 


81 


PKS2155-304 


1238.56 


0.22 


52 


28 


39 


36 


33 


PKS2155-304 


1222.59 


0.32 


78 


384 


241 


56 


42 


Q1230+0115 


1222.98 


0.46 


113 


241 


141 


42 


50 


Q1230+0115 
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7. Metallicity of the Low-z Lya Forest 

One of the interesting measurements that may elucidate the high-z origin of the low-reelshift 
Lya absorbers is a search for heavy elements. Weak C IV lines are detected at high-z in 50-75% of 
the Lya forest clouds with Ahi > 10 14 5 cm~ 2 (Songaila & Cowie 1996) and perhaps to even lower 
column densities (Ellison et al. 2000). The C IV/H I ratios suggest abundances of io _2 - 5±0 - 5 times 
solar, and the Si IV/C IV ratios suggest that these elements have been formed by early generations 
of massive stars and expelled by galactic winds or tidal stripping (Giroux & Shull 1997; Gnedin 
1998). Observing the metallicity of the low-z absorbers would help to understand the chemical 
evolution of the IGM, and would test the prediction from numerical simulations (Cen & Ostriker 
1999b) that the IGM metallicity has increased to values of 10-20% solar at the current epoch. 

Our first metallicity observations were performed on the cluster of strong Lya absorbers toward 
PKS 2155-304 (Shull et al. 1998). ORFEUS observations along this sightline (Appenzeller et al. 
1995) suggested, from Lyman continuum absorption, that A^hi = (2 — 5) x 10 16 cm~ 2 at the redshift 
of the Lya lines. The saturated Lya lines only provide a conservative range of (3 — 10) x 10 14 cm~ 2 
from our HST/GHRS data. Our recent measurements of higher Lyman series lines and the Lyman 
continuum edge with FUSE show that Nm < 10 16 cm~ 2 . From the non-detection of Si III A1206.5 
and C IV A1548.2 at their predicted positions, together with photoionization corrections, we hope 
to improve the limits on (Si/H) and (C/H) using HST/STIS spectra obtained in cycle 8. 

To explore the mean metallicity of the low-z Lya forest, we constructed a composite spectrum 
for all Lya absorbers. Using a "shift and stack" method, we aligned all SL > 4a Lya absorbers 
in the absorber rest frame and then co-added the spectra. Virtually all of these Lya absorbers lie 
well away from any predicted metal-line position. However, we have not co-added any portion of 
the spectrum whose wavelength region (5 A around the metal-line position) included a definite Lya 
absorber. In one case, the Lya absorber at 1226.964 A (2785 km s" 1 ) towards PKS 2155-304, 
has a slight chance of being Si III A1206.500. Although we believe this alternate identification is 
unlikely, based on photoionization and metallicity considerations, we have not included it in our 
coaddition. 

Table 7 indicates the strong metal resonance lines in the wavelength region of our spectra and 
the upper limits on W me t a is for each species. All reported W values are 4<r detection limits in mA, 
based upon the pixel-to-pixel variations in the region of the expected metal line of the rest-frame 
composite absorber spectrum. Because the composite spectra include our pre-COSTAR observa- 
tions, the determined 6- values of Lya are slight overestimates due to the spectral smearing of these 
observations from target motion in the GHRS LSA (§ 2). We divide our definite Lya sample into 
two sub-samples based upon the equivalent width (Wiya) of the Lya absorber. Table 7 gives metal- 
licity limits, in the low W sample (Wl ?/q ,<133 mA), in the stronger absorbers (Wl ?/q ,>133 mA), 
and in the combined sample (column denoted "all" ) . This division is designed to separate saturated 
and unsaturated absorption lines based upon their n(W) distribution (Figure 7). Searching the 
available rest-frame spectral regions (1138 A < A < 1298 A) for possible metal lines, we found no 
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absorption down to the 4a level in any W ' Lya sub-sample. 

The composite ultraviolet spectrum can be used to place an upper limit on the average metal- 
licity associated with our set of Lya absorbers. Unfortunately, the strong carbon lines, C III A977 
and C IV A1549, lie outside our A-range. Most of the resonance lines of Table 7 have ionization 
thresholds similar to that of hydrogen, except for Si III A1206.5 which is 33.5 eV. The lines with 
~15 eV ionization thresholds have ionization corrections (IC) similar to that for hydrogen, which 
are quite high, resulting in low predicted Ws. Si III has a much lower IC, which equates to more 
Si III available for detection and a higher predicted W(A1206.5). Therefore, of the resonance lines 
in our waveband (see Table 7), our most promising limit on the metallicity is that derived from 
Si III, due to its higher ionization threshold. 

Table 7: Metallicity upper limits (4cr ) for our composite low-z Lya absorber 



Feature 


A(A) 


All 


W metals (Definite Sample) 
W Lya < 133mA W Lya > 133 mA 


SiHI 


1206.50 


2 


2 5 


Mnll 


1201.12 


2 


2 7 


NI 


1200.71 


2 


2 7 


NI 


1200.22 


2 


2 7 


NI 


1199.55 


2 


2 8 


Mnll 


1199.39 


2 


2 8 


Sill 


1197.39 


2 


3 12 


Mnll 


1197.18 


2 


3 12 


Sill 


1193.29 


3 


4 18 


Sill 


1190.42 


4 


5 22 



The composite Lya spectrum for our low-W, unsaturated sample (Wz^qX 133 mA) is shown 
in the upper left panel of Figure 22. The measured W is 41±2 mA and the measured 6-value of 
the composite absorption is 6 b s = 38±1 km s" 1 . No Si III A1206.5 absorption is detected in the 
low-W composite spectrum down to the 2 mA (4a) level. The composite Si III A1206.5 spectrum for 
the low-W sample is shown in the upper right panel of Figure 22. In our sample of unsaturated Lya 
forest lines, the ratio of the composite equivalent widths is W(Si III A1206.5)/W_L ?/a < 2 mA/41 mA 
~0.05. For unsaturated lines, the composite absorption should be a good approximation of the 
average absorber since W scales linearly with Nni- 

In the linear regime, the rest-frame equivalent width (W) is related to the column density of 
absorbing atoms along the line of sight, Nj, by: 

2 

W (cm) = ^-X 2 Njf jk = 8.85 x l(T 13 A 2 A^ fc , (22) 

where Nj is in cm" 2 , fj^ is the oscillator strength of the transition, and A is the wavelength of the 
absorbed photon in cm (Spitzer 1978). With current oscillator strengths (0.4162 for Lya; 1.68 for 
Si III A1206.5, Morton 1991) and wavelengths, the relationships between W and AT HI for these 
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transitions are: 

W Lya = 54.43 



10 13 



cm 



mA, and Wgiin = 2.164 ( J^ S ™_ 2 ) mA , (23) 



which translates to a ratio of column densities of: 



iVsnii Q2514 fWsim , <0012) 



N m v whi ; ' (24) 

given the observed limits. Determining the actual cloud metallicities (Z) requires knowledge of 
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Fig. 22. — Composite Lya and Si III A1206.5 spectra for all absorbers with W Lya<133 mA (upper 
panels) and V^Lya> 133 mA (bottom panels). 

the fractional amounts of Si III (/ci = Si Ill/Si) and H I (/hi= H I/H) in the absorbing clouds 
(Donahue & Shull 1991; Shull et al. 1998). It is customary to express metallicities (Z) relative to 
solar (0) values in logarithmic form (denoted by [ ]). The metallicity of the absorbing clouds is 



Z Q 



Si 
H 



N S i 



/hi 



VhI ^ /Silll ' 



(h)g (h% 



or 



(25) 
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"Si 
H 

The ionization correction (IC) in this form accounts for the ratio of fractional abundances corrected 
for the solar abundance ratio. For Si, the solar abundance ratio, log(Si/H) Q = —4.45 (Grevesse & 
Anders 1989). The ionization correction (IC) can be made and the metallicity of the clouds deter- 
mined, only if /siin and /hi are known. Applying the solar abundance of Si and our observational 
limits equation (24) to equation (27) provides: 

< log (0.012) +IC = -1.91 + 4.45 + log ( ) = 2.54 + log ( J . (28) 

V/SiIII/ V/SiIII/ 

The determination of the fractional abundances depends on uncertain photoionization modeling of 
the gas responsible for the absorption. We have used the photoionization code CLOUDY (Ferland 
1996) to calculate the fractions of H I and Si III A1206.5 expected for three representative models 
of the gas. The measured W and 6-value of the low-W composite spectrum indicate a composite 
A%i = 8x 10 12 cm -2 . All models assume that the gas is optically thin, illuminated by a metagalactic 
radiation field with a power-law spectrum with spectral index a s = 1.8, and mean ionizing intensity 
Jo = 10~ 23 ergs cm" 2 s" 1 Hz -1 sr" 1 . We consider three models with total hydrogen densities 
nn= 10~ 4 ,10" 5 , and 10~ 6 cm" 3 . The resultant ratios of log (/m/Zsiin) are —2.4, —1.6, and +1.5, 
respectively, implying [Si/H] <0.14, 0.94, and 4.02 solar, or Z Si /Z Q < 1.4, 8.7, and 10 4 . The very 
high upper limit on metallicity for the njj = 10~ 6 cm" 3 model is due to the transference of the 
dominant state of Si to Si IV. The characteristic density for these absorbers remains an unknown 
quantity, and it may be seen from this calculation that a strong constraint on metallicity is not 
yet possible regardless of the cloud density. In Shull et al. (1998) we were able to place restrictions 
upon i%h along the PKS 2155-304 sightline by assuming cloud depths that are constrained by VLA 
observations of nearby H I galaxies. Typically, observations of two or three separate states of a 
single element (e.g., C II, C III, and C IV) are required to measure the metallicity accurately. 
Therefore, the limited wavelength coverage of our observations prevented a more stringent limit of 
Lya cloud metallicity. 

The composite spectrum for our high-W, saturated, Lya sample (Wlj /q >133 mA) is shown 
as the lower left panel of Figure 22. The measured W is 231±4 mA and the measured 6-value of 
the composite absorption, adjusted for the resolution of the GHRS/G60M, is 6 b s = 51±1 km s -1 . 
As shown in the lower right panel of Figure 22, no Si III A1206.5 absorption is detected in the 
Wz /2/a >133 mA composite spectrum down to the 5 mA (4a) level. Because some absorbers at 
W > 133 mA may be saturated, the composite value of A%i is uncertain. If we adopt 6=25 km s" 1 , 
then the fiducial W (231 mA) yields N m = 7.6 x 10 13 cm , whereas we obtain Nm = 3.9 x 10 
cm" 2 for an extension of the linear curve of growth (equation 23). For Si III A1206.5, as determined 
from equation 23, Nsmi <2.3 x 10 11 cm" 2 . Applying the same ionization correction as for the low- 
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W composite spectrum, we estimate: 



[I] < * (SH = + « 5 + * (&) — + * (&) • <»> 

For the three CLOUDY models, with assumed hydrogen densities of nn= 10 _4 ,10 -5 , and 1CT 6 
cm" 3 , we find [Si/H] < -0.47,0.33, and 3.43 solar, or Z Si < 0.33,2.13, and 3 x 10 3 Z Q . Again, 
these limits are not at all stringent because absorption lines from only a single ionization state of 
Si fall within our observed wavebands. 



8. The Contribution of the Local Lya Forest (FlLya) to fl^. 

The contribution of the low- z Lya forest to the total number of baryons, as estimated from 
Big Bang nucleosynthesis, can be determined from: 

a Lm = f" " , ( 30) 

J N mm Per 

where </>o(-/Vhi>p) is the comoving space density of the clouds, M c i(AT HI ,p, J ) is the individual cloud 
mass, p is the impact parameter of the line of sight from the cloud center, TVhi is the observed column 
density, p cr = 3/Jq / 8nG is the present-day critical density necessary to halt the expansion of the 
Universe, Jo is the specific intensity of the metagalactic ionizing radiation field at hvo = lRyd, 
in ergs cm -2 s _1 Hz -1 sr _1 , and (N min ,N max ) is the range of integration appropriate for the 
assumed cloud model. Assuming that all the low-z Lya absorbers are singular isothermal spheres 
with density distributions of total hydrogen, nn(r) = no(r/ro)~ 2 , the cloud mass in the absence of 
metals, but including helium, within radius r = p is given by: 

rp 

M cl (N m ,p, J ) = / n H (r)nb^TTr 2 )dr = A7rn Q rl(l.22m H )p , (31) 
Jo 

where \x\, = 1.22m# is the average mass per baryon, assuming that the fraction of helium by mass 
is Y = 0.24 or x = n He /n H = 0.0789. 

Assuming that the gas is in photoionization equilibrium at T = 20,000 K, immersed in an 
ionizing radiation field, J u = Jo (u/vo)~ aa , we can equate the observed column density to the 
expected H I column density by integrating through the cloud for a given impact parameter, p, 

/•oo 

N m (p) =2 /hi n H (r) dl . (32) 
Jo 

Here / is the distance along the line of sight measured from the point of closest approach to the 
cloud center, p, so that I 2 = r 2 — p 2 . The ionization fraction of H I is given by: 



-52- 



where a^p , the case-A hydrogen recombination rate coefficient, 3 is 2.51 x 10 13 cm 3 s 1 at T = 
20, 000 K, and the photoionization rate, Thi, in the presence of the described radiation field, J u , is 

= 4, (£) = tS^T - (2-49 * 10- C " 1 , (34) 







hv J h(a s + 3) \a s + 3 



where a s ~ 1.8 (Zheng et al. 1997), J-23= Jo / (10 ergs cm s Hz sr~ ), and a v m 
do {v/vq) 3 is the hydrogen cross section at frequency v. The scaling of Jo is supported by a recent 
theoretical estimate of Jo = 1.3^q g x 10~ 23 ergs cm~ 2 s _1 Hz -1 sr _1 in the local Universe (Shull 
et al. 1999b) and the upper limit of Donahue, Aldering, & Stocke (1995) of J < 4.2 x 10~ 23 ergs 
cm" 2 s _1 Hz -1 sr _1 for a s = 1.8. Integrating through the cloud, using nm = n e nH+a^f V^m and 
n e = (1 + 2%) gives the column density, 

/°° In cv^ f°° 

n ffl (r) dl = (1 + 2 X ) / [n ff (r)f d/ (35) 

-oo 1 HI Jo 

vrn 2 r 4 «^ ) (l + 2 X ) 

2r ffl p 3 • (36) 

Spectroscopic studies of multiple QSO sightlines at high-z (Bechtold et al. 1994), and nearest- 
neighbor studies of low-z Lya absorbers (Stocke et al. 1995; Shull, Stocke, & Penton 1996; Shull 
et al. 1998) indicate that Lya absorbers may have spatial extents of ~100 kpc. However, the sight- 
line pair experiments have found significantly larger values for individual Lya absorbers (Dinwhaw 
et al. 1997). We adopt this as a conservative impact parameter so that pioo=p/(100 kpc), recog- 
nizing that some of the Iow-Ahi absorbers in voids may be considerable larger. Solving equation 
(36) for norg, and combining this result with equation (31), gives the Lya cloud mass as : 



M cl (N m ,p, J , a s ) =47r(1.22m JT )p5 



2r H iA^Hi 

(1 + 2 X ) 



(37) 



(1.59 x 1O 9 M ) Nl J^ 23 4o (^fs) 2 > ( 38 ) 



where Ni4=Nm/ (10 14 cm 2 ). At each iVni, the space density of absorbers, 4>q, is related to oW /dz 
by: 

H <W{N m ) 
crrp 2 dz 

Therefore, assuming that p(Nm) is a constant 4 , we find 



MNm) = 7T^2 JZ^ 1 ■ (39) 



ya 



(1.59 x 1O 9 M ) H J^ 23 pj 00 (^) 2 f N maX(W{Nm) 1 

2 / j iV 14 dAr H i , (40) 

cnp p cr J Nmm dz 



3 Case-A recombination assumes that the energetic photon emitted in ground-state recombination is not reabsorbed 
by the Lya cloud. 

4 The low- z clouds in galaxy voids could have systematically larger impact parameters for a given column density, 
making the impact parameter versus column density relationship more complicated. 
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where, for an flo = 1 Universe, p cr = (9.21 x 10~ 30 ) /i 2 g cm -3 . Performing the integration 
of equation (40) for the low-z Lya forest over the high end of our column density range, 13.5 < 
log [Nffl] < 16.0 gives: 

• (41) 

This value represents ~20% of the baryons (Q^hj = 0.0388) as inferred from D/H by Buries & 
Tytler (1998). Decreasing N min to 10 123 cm -2 would increase this estimate by roughly a factor of 
2 to £lLya = 0.018 ± 0.003. However, the absorbers with log [A%i]<13.5 are probably uncollapsed 
matter, for which isothermal spheres of constant impact parameter are increasingly poor physical 
approximations. In addition, the assumed power-law distribution in column density may not extend 
to such low values of Ahi- 

There are various corrections to our simple model that could be applied to our estimates. 

Numerical simulations (Dave et al. 1999) suggest that a flattened geometry is more appropriate for 

the higher column density Lya systems (Dinwhaw et al. 1997; Crotts & Fang 1998). Thus, the 

absorbers may be denser than assumed, with smaller ionization corrections. As pointed out in Rauch 

& Haehnelt (1995), Shull, Stocke, & Penton (1996), and Madau & Shull (1996), by generalizing 

i i_ 

from spheres to disk geometries, the value for £liya is reduced by a factor of (acos#)~2 ~ (2/a) 2 , 
where a (>1) is the disk aspect ratio and 9 is the disk viewing angle ((cos 6) = 1/2). For example, 
an aspect ratio of 10:1 would reduce VLiya by a factor of 2.2. Further complicating the conversion 
from Ahi to Q,Lya is the possibility that a substantial portion of the low-z Lya clouds are hot 
(10 5 — 10 7 K) and collisionally ionized (Cen & Ostriker 1999a). Thus, our ionization correction 
from H I to total H may be too low for some clouds, and others may be undetectable given our 
sensitivity limits. Although the results of the ^Lya contribution to fib are model dependent and 
highly uncertain, it is likely that warm, photoionized Lya clouds contain a substantial fraction of 
the baryons, &Lya= 0.008-0.018, or 20-45% of Q^. This agrees with numerical simulations that 
suggest that a substantial fraction of baryons in the local Universe lie in these Lya absorbers. 



n Lya = (0.008 ± 0.001) 



J-23 P100 



4.8 



o , 



+ 3 
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9. Conclusions 

The major conclusions of our physical analysis for the low- 2; Lyo forest are: 

• Our observed 6-value distribution is consistent with that at higher redshift. We see little 
evidence for an increasing median 6-value with decreasing z as reported by Kim et al. (1997). 
We are cautious about inferring that our measured 6-values (from line widths) are accurate 
indicators of the true Doppler parameters of the low- z Lyo forest. The unknown contribution 
of intracloud turbulence or the presence of multiple sub-components per absorber leads us 
to adopt a constant value 6=25 km s _1 when converting equivalent widths (W) into column 
densities (Nm)- 

• Our observed W distribution appears similar to that in higher redshift Keck/HIRES data. 
Applying a careful sensitivity correction, we detect a significant break in oW /dW at W < 
133 mA, which probably arises from saturation in the Lyo curve of growth. Our results for 
dN /dz are consistent with HST/FOS Key Project low-z studies, but they extend to much 
lower equivalent widths, W~14 mA for 4<r and W~10 mA for 3<r detections. 

• Using an integrated method (Table 5) for evaluating the functional form of d 2 Af/dz dNm, 
we find that, for log [N m ] < 14, d 2 M/dz dN m oc N u f with (3 = 1.81 ± 0.05. There is some 
evidence for a break at log [Nm] = 14, above which (3 = 1.43 ± 0.35. The reality of this break 
is suspect due to small number statistics and to the dependence of column density on our 
selection of 6 = 25 km s _1 . 

• We find no z evolution over our small redshift range (0.002 < z < 0.069), and we observe 
no unusual deviations in dAf /dz that could arise from large-scale structure in the local 
Universe. We calculate that the low-z Lyo forest produces a Lyman-continuum opacity at 1 
Ryd of dr cff /dzRi0.01 for log [N m ] < 13 and ~0.1 for log [N m ] < 15. Most of the intergalactic 
opacity probably arises from rare, higher- 7Vhi absorbers. 

• We have compared the z evolution of oW /dz (Figure 17) and d 2 N /dz cWhi (Figure 18) at 
low and high iVni and z. In both cases, we find significant evolution that can be explained 
by: (1) expansion of the Universe, (2) the rapid decline in J v {z) at z<2, and (3) the merger 
and dissipation of Lyo absorbers. Improved statistics over intermediate redshifts (0.1 < z < 
1) and a wide range in Nm are required to extract more meaningful evolutionary results. 
However, these preliminary results agree with hydro dynamical simulations of the low-z Lya 
forest by Dave et al. (1999). 

• By comparing the z evolution of dj\f /dz for the low column density (13.1 < log [Nm] < 
14.0) absorbers studied in this paper with similar results for higher column density absorbers 
(log [Nm] > 14) found by the Key Project team, we find no difference in the overall z evolution 
of the low and high Ahi absorbers. This is in disagreement with previous results at both 
high and low-z (Kim et al. (1997) and Weymann et al. (1998), respectively) who found slower 
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evolution for the low column density absorbers. However, lacking information at intermediate 
redshifts, all these results can be in agreement if the dramatic "break" in dA/" /dz occurs later 
(i.e., at lower redshift) for the lower column density absorbers. 

• We find a 4cr signal in the two-point correlation function (TPCF) of \ow-z Lya absorbers for 
velocity separations Av < 250 km s _1 , consistent with results from higher redshift studies. 
Our results at higher velocity separations are limited by relatively poor statistics. In Paper IV, 
we will combine these GHRS observations with our results from 13 STIS sightlines to improve 
our understanding of the low-z Lya TPCF. 

• Applying a photoionization correction, we find that the low-z Lya forest may contain ~ 20% 
of the total number of baryons, with closure parameter Qi,ya = (0.008 ± O.OOl)/^ 1 , for a 
standard absorber size of 100 kpc and an ionizing radiation field of intensity Jo = 10~ 23 ergs 
cm~ 2 s _1 Hz -1 sr _1 and spectral index a s = 1.8. 
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staff at the Space Telescope Science Institute, particularly Ray Lucas. We thank Mark Giroux for 
helpful discussions and a critical reading of the manuscript. Buell Jannuzi and Ray Weymann are 
thanked for a critical reading of the manuscript and helpful conversations. This work was supported 
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Astrophysical Theory Program (NASA grant NAGW-766 and NSF grant AST96-17073). 
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